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Abstract. The following review is based on lectures given in the 1 st Samos Cosmology 
summer school. It presents an attempt to discuss various issues of current interest 
in Observational Cosmology, the selection of which as well as the emphasis given, 
reflects my own preference and biases. After presenting some Cosmological basics, for 
which I was aided by excellent text-books, I emphasize on attempts to determine 
some of the important cosmological parameters; the Hubble constant, the curvature 
and total mass content of the Universe. The outcome of these very recent studies 
is that the concordance model, that fits the majority of observations, is that with 
Q m + Q A = 1, i^A ~ 0.7, H ~ 70 km s" 1 Mpc -1 , J? B ^ 0.04 and spectral index of 
primordial fluctuations, the inflationary value n ~ 1. I apologise before hand for the 
many important works that I have omitted and for the possible misunderstanding of 
those presented. 

1 Background- Prerequisites 

The main task of Observational Cosmology is to identify which of the idealized 
models, that theoretical Cosmologists construct, relates to the Universe we live 
in. One may think that since we cannot perform experiments and study, in a 
laboratory sense, the Universe as a whole, this is a futile task. Nature however 
has been graceful, and through the detailed and exhaustive analysis of the de- 
tected electromagnetic radiation emitted from the different photon-generating 
processes, we can do wonders! 

Among the many important tasks of Observational Cosmology is the determi- 
nation of the total mass-energy density of the Universe, the rate of its expansion, 
its age, the amount of ordinary and exotic matter that it contains, as well as 
to quantify in a objective and bias free manner the large-scale distribution of 
matter, as traced by galaxies, clusters of galaxies and AGN's for example. 

However, these tasks are not easy to fulfil. Subjective (instrumentation, avail- 
able funds, technological limitations etc) as well as objective (observational bi- 
ases, limitations due to our position in space-time, etc) difficulties exist. Further- 
more, we do not know whether the Universe accessible to our observations is rep- 
resentative of the whole Universe. A positive answer to this question is essential 
in order to meaningfully compare observations with theory. Under the assump- 
tion that the Universe is homogeneous and isotropic (in a statistical sense), well 
separated regions can be viewed as independent realizations of the same forma- 
tion process. Therefore many of such regions constitute an ensemble and thus 
we can employ statistical techniques in our study. 
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In the classical Big-Bang cosmological framework the Universe indeed is 
considered homogeneous and isotropic on the large-scales. The most general 
metric satisfying this assumption, the so-called Cosmological Principle, is the 



Robertson-Walker metric (cf. [183|, |35 ), [112Q 



ds 2 = c 2 dt 2 - R 2 {t) 



dr 



r 2 {d6 2 +sm 2 ed<j) 2 ) 



(1) 



where R{t) is the expansion factor, A; is a constant, related to curvature of space 
and (r, 9, <f>) are spherical-polar coordinates. The main observational evidence 
that supports the choice of this model is: 

• The observed expansion of the Universe. Edwin Hubble found that the red- 
shifts of galaxies are proportional to their apparent magnitudes and assuming 
that they are equally luminous then their redshifts are proportional to their 
distances: v cx d. 

• The existence of the cosmic microwave background (CMB) radiation, inter- 
preted as the relic radiation from the hot initial phase of the Universe. 

• The observed light element abundances that this theory correctly predicts. 

However, the observed matter distribution in the Universe is very inhomogeneous 
on small-scales. So, what evidence do we have supporting the validity of the 
Cosmological Principle? 

Firstly, the Hubble's law is directly obtained if one assumes a homogeneous 
expansion of the Universe, ie., if a length \ is expanded by a factor R(t), then 
after some time we have d — R{t)\- Differentiating we obtain the Hubble law: 

v = ^d = H(t)d, (2) 

where H(t), at the present time (t = 0), is the Hubble constant. Secondly, 
observations of distant extragalactic radio sources have shown that they are dis- 
tributed across the sky in a uniform way. Other supporting evidence is provided 
from the decreasing - with scale - correlations of extragalactic objects and the 
directional independence of the correlation function. However, the most remark- 
able confirmation that the Universe is homogeneous and isotropic and also that 
it evolved from a hot dense past, was the discovery of the cosmic microwave 
background radiation and its high degree of isotropy across the sky. This ra- 
diation has been interpreted as the relic radiation from the time that matter 
decoupled from radiation, which has been freely travelling ever since. The high 
degree of isotropy of the CMB is direct evidence that the Universe was highly 
isotropic at the epoch of decoupling (z ~ 1100) to one part in 10 5 on scales from 
arc-minutes to 90°, once we subtract a local dipole anisotropy, attributed to our 
peculiar motion with respect to the rest-frame defined by the CMB. 



1.1 Basic Elements of Dynamical Cosmology 

Out of the many parameters which are essential for accurately determining the 
global dynamics and the formation history of the Universe; three clearly stand 
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out: the Hubble constant, H , which tells us the expansion rate of the Universe, 
the matter density parameter, J? m , which tells us how much matter, being bary- 
onic or exotic, the Universe contains and the Cosmological constant, A, which 
tells us whether the universe is filled with an extra repulsive force. The way in 
which these parameters affect the cosmological evolution, are determined by the 
gravitational field equations and the assumed equation of state. Below, I will 
sketch the basic framework of standard Cosmology, in order to derive the inter- 
relations of these parameters and how they relate to the global dynamics of the 
Universe. 



Priedmann Equations: Within the context of a homogeneous and isotropic 
Universe we can derive the cosmological evolution equations using either Ein- 
stein's field equations or Newtonian gravity. The latter is possible exactly be- 
cause of the Cosmological Principle, ie., we can consider any volume element, as 
small as necessary for Newtonian gravity to apply, as being representative of the 
Universe as a whole. We will derive the evolution equations for the Newtonian 
case, using however as an active mass density: 



P = p + ip/c 2 



(3) 



which is given within the framework of general relativity for a homogeneous and 
isotropic fluid with density p and pressure p. 

In this case we can derive the evolution equations using the mass continuity 
equation and Newton's equation of the motion of a small sphere. We have from 
the homogeneity assumption that V p = 0, from isotropy that V • v = 3H(t) = 
3R/R and then from dp/dt + V • (pv) = we obtain: 



p + 3 



P 



R 
R 







Rearranging (Q) and using Newton's equation of motion, R 
have: 



d_ 

df 



1 



:R 2 



4nG 



pR 2 



= 



(4) 



-GM/R 2 , we 



(5) 



and integrating we obtain a Newtonian analogue of cosmological evolution equa- 
tion: 

R 2 8ttG C 

(6) 



R 2 



where C is the constant of integration, which is closely related to the Newtonian 
energy. We can see this if we rearrange @ as: R 2 /2 — GM/R = —C/2. Within 
the relativistic formulation, the constant is in effect related to the curvature of 
space; C = kc 2 and the basic equation of cosmological evolution, the Friedmann 
equation, is written: 

R 2 kc 2 _ 8irGp Ac 2 
R 2 + "R 5 " - 3 + — 



(7) 
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where the cosmological constant (A) term was introduced by Einstein (rather ad 
hoc) in order to obtain his preferred static solution. Given an equation of state, 
P = p{p)i we can solve for p and then (Q) can be integrated to give R(t). 

If we recall Newton's equation of motion, R — —GM/R 2 , and we use the 
active density (|^) and (Q), we derive the second important dynamical equation 
of cosmological evolution. The correct relativistic form of this equation is: 

2R R 2 kc 2 8nGp , , 

-R + R 2+ 1?=-— +AC < 8 > 



Equation of State: The question arises of which is the appropriate equation 
of state for the expanding Universe. The Universe expands adiabaticallyQ, since 
the symmetry imposed by the Cosmological principle implies that there is no net 
heat flow through any surface. Therefore as it expands, it cools ( |Ti"| ) and since it 
started with a very hot phase, it will be dominated, at different epochs, by differ- 
ent species of particles having distinct equations of state. Relativistic particles 
will dominate the hot phase while non-relativistic the later cooler phases. 

We can specify a unique equation of state p — p(p) for all epochs by param- 
eterizing it according to: 

p = w(v 2 )p (9) 

where (v 2 ) is the velocity dispersion of the fluid elements. If the dominant con- 
tribution to the density comes from relativistic particles, which havep = l/3pc 2 
(when kT 3> m v 2 , with v 2 ~ c 2 ), then w — 1/3. If the dominant contribution 
comes from non-relativistic matter (when kT ~ m v 2 with v 2 <C c 2 ) then there 
is negligible pressure and the dust approximation is excellent (w = 0). 

Therefore inserting (^J) into the mass continuity equation, (^), we obtain: 

p cx ir 3(1+tu) (io) 



Armed with (10) one can now solve the Friedmann equation to get the time 
evolution of R(t), determine the age of the Universe, etc, in the different Cos- 
mological models. 



1.2 Thermal Beginning of the Universe 

The early universe, where very high densities and temperatures dominate, can 
be treated by using fluid thermodynamics. At very high temperatures, radiation 
and matter are in thermal equilibrium, coupled via Thomson scattering with the 
photons dominating over the nucleons (n 7 /n p ~ 10 9 ). Therefore the primordial 
fluid can be treated as radiation-dominated with p — l/ipc 2 = 1 /3(jT 4 and from 
([lO|), we obtain: 

TtxR- 1 (11) 

1 An adiabatic process is defined as that in which there is no heat flow and thus the 
entropy is conserved (dS = 0) . 
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Therefore the temperature of the Universe drops linearly with the expansion 
scale factor. Furthermore, it is evident from (flO|), that the radiation density 
drops faster than the mass density and since we know from measurements that 
the universe is matter dominated today, then at some epoch in the past, say 
at a redshift z cq , we had p m = p ra d- It is easy to show that p T = p m R /R cq = 
(1 + z cq ) p nl (the subscript o denotes the present epoch) and using the measured 
values of pi we have that: 

1 + z cq ~ 2.3 x 10 4 h 2 Q m 

Therefore the thermal history of the Universe can be divided in two main eras: 
a radiation dominated era (z ^> z cq ) and a matter dominated era (z <§C z cq ). In 
the radiation dominated era, in which we can neglect the curvature and A terms 
in Fricdmann's equation (see next section), we have: 

Rozt 1 ' 2 . 

By differentiating this relation with respect to time and using (^) we have: 

1/2 

(12) 



327rGp 7 

Using p~ f = TT 2 kbT 4 /15h 3 c 5 we finally obtain the important relation between 
cosmic time and the temperature of the Universe in the radiation dominated 
era: 

T Kclvin ~ 1.3 x I0 10 t-V 2 (13) 

from which it is evident that the Universe at early times was hot enough for 
nucleosynthesis to occur, as it had been supposed originally by Gamow. The era 
of nucleosynthesis takes place around ~ 10 9 K. 



The Cosmic Microwave Background: Although the dynamics during the 
radiation dominated era are unaffected by ordinary matter, the electrons act 
as a scattering medium of the radiation and thus the Universe at this epoch 
is opaque. As the Universe cools, oc electrons bind electrostatically with 

protons to form Neutral Hydrogen. Using Saha's ionization equation one finds 
that the temperature at which the number of free electrons drops significantly 
is T ~ 3000 K. 

Therefore when the universe cools at this temperature, the scattering medium 
disappears and the radiation freely escapes without being absorbed or scattered 
which means that the Universe becomes transparent. This epoch is called the 
recombination epoch. 

The existence of the relics of this radiation was predicted by Gamow and his 
collaborators in the 1940s. It was subsequently discovered by Penzias & Wilson 
in 1965, while the whole spectrum of this radiation was traced to unprecedented 
accuracy by the COBE satellite observations. The CMB possesses a perfect 
black-body spectrum with a mean temperature of T — 2.728 ± 0.004 K and it 
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is extremely isotropic except for a dipole, which is however a local kinematical 
effect (due to our motion with respect to the cosmic rest frame defined by the 
CMB). From what redshift does the CMB radiation originate? From ( pd| ) we 
have that: 

T Ro 

T =ll =1 + Z 
with T ~ 3000 K and T ~ 2.73 we get that 

l + z rcc ~ 1100 

From (0) and (|To|) we have that in the matter dominated era R(t) oc ti and thus 
z rcc corresponds to a time: 

i rcc ~ 2.8 x 10~ 5 t 

where t Q is the present age of the Universe. Therefore by studying the microwave 
background sky we have direct information from the Universe when it was as 
young as t Tec . 



The CMB dipole anisotropy: Due to our motion with respect to the isotropic 
CMB radiation we observe a dipole in the distribution of the radiation tempera- 
ture. Although this has the appearance of a Doppler effect, in reality four different 
effects add up to produce this dipole seen by an observer moving with a velocity 
u. These four effects are: 

• a Doppler effect that increases the frequency of photons, and thus the ob- 
served energy, seen in the direction of motion by a Doppler factor D = 
1 + (u/c) cos 9 

• the interval of frequencies increases by the same factor in the direction of 
motion, and therefore since T oc E/Sv, the above two effects cancel out. 

• the moving observer selects in the direction of motion relatively more photons 
by a factor D 

• the solid angle in the direction of motion is smaller by a factor D~ 2 due to 
abberation. 

The net effect is that the moving observer sees an intensity of CMB radiation 
Jmov = (1 + u/c cos 9) 3 I lcat . Due to the adiabatic expansion of the Universe, 
(T oc R^ 1 ), the shape of the Planck spectrum: 

m-r 

should be preserved, which then necessarily implies that T(9) = (1+u/c cos 9)To 
and thus: 

AT u 



h = 
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COBE observed a CMB dipole amplitude of ST ~ 3.3(±0.2) mK (which corre- 
sponds to a fluctuation ST/T = 1.2(±0.03) x 10~ 3 ). The corresponding velocity 
of Earth is: 

Vq — vcmb ~ 365(±18) km/sec 



towards the galactic coordinates (l,b) = (265°, 48°) (see |159|| ). This motion is 
due to the vectorial sum of the motion of the Earth around the Sun, of the Sun 
within the Galaxy, of the Galaxy within the Local Group and of the peculiar 
motion of the Local Group, due to the gravitational effects of large-scale density 
fluctuations. The motion of the Earth with respect to the LG centroid is: 

Vq — i>lg ~ 308 km/sec 

towards (l,b) = (107°,— 7°) and thus we find the velocity of the LG centroid 
with respect to the CMB: 

vlg - v C mb ~ 620 km/sec 

towards (l,b) = (277°, 30°) . 

The Local Group velocity was originally thought as the result of the attrac- 
tion of the Local Supercluster (Virgo). However, there is a residual velocity of 
~ 400 km/sec that must be due to gravitational forces acting on the LG from dis- 
tances greater than the Local Supercluster's centre-of-mass (cz ~ 1100 km/sec). 
Many earlier studies pointed towards the 'Great Attractor', a mass concentra- 
tion of ~ 5 x 10 16 Mq located at a distance of 42 h~ x Mpc and at low Galactic 
latitudes, as being the sole cause of a relatively local coherent motion, in which 
the Local Group partakes (cf. [pi| , [p3|). Later studies, indicated that another 
very massive and more distant (~ 140 h^ 1 Mpc) attractor could play a signif- 



icant role in shaping the local dynamics ([151], [152], ]122|] ). It seems that the 
coherence scale of the velocity field could extend to even larger distances than 
what originally thought (cf. O], however for a different view see pOj). 



1.3 Cosmological Parameters 

Based on (0) we can define some very important parameters like the Critical 
Density, which is the density necessary to obtain a flat Universe (A = 0): 

Per = ¥^ = L88 x l(T 29 /i 2 gm cm" 3 (15) 
87rG 

and the Cosmological density parameter fl, which is a unit-less measure of the 
density of the Universe: 

n = JL (16) 

Per 

Furthermore, the constant of proportionality in Hubble's law, the Hubble con- 
stant, is: 

Ho = 100 h km = 1.023 x 10~ 10 h years^ 1 (17) 
sec Mpc J v ' 
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Note that the necessity of parametrizing with h was due to earlier discordant 
determinations of H . Today most studies converge to a value of ~ 0.7 (see 
section 2.4). 

A convenient representation of these interrelations can be produced by re- 
writing Friedmann's equation, in the matter dominated era (using |l0| ), as fol- 
lowing: 

| = H (n m (i + zf + n k (i + z f + n A ) 1/2 h{z) = h e(z) (is) 

were the contribution to the total density parameter from the curvature and A 
terms is: 

kc 2 Ac 2 

Ok—HM' (19) 

Note that H(z) is called Hubble function. It is evident that at the present epoch 
we obtain from ( |l8| ) that E(0) = 1 and thus: 

Q m + f2 k + f2 A = 1 (20) 

which also holds for any epoch (evaluated directly from 0). Note that we can 
have a flat Universe (J7k = 0) while having f2 m < 1 (as suggested by many 
different observations). 

The Age of the Universe: Using (|l^), evaluated at the present epoch, we 
have R/R = H E(z)/(l + z) and from dR/R = —dz/(l + z) 2 we obtain the 
age of the Universe: 

to -E-J Q (l + z)E(z) (21) 
For example, in an Einstein-de Sitter universe (f2\ — f2 k = 0) we have: 

f ° = dk (22) 

while for a J?a > model we obtain: 

91 . r nrr 

(23) 



t^ = — -Lsinh" 1 

2>H y 



We therefore see that if ]?a > we have that the age of the Universe is larger 
than what is predicted in an Einstein-de Sitter Universe. 



The A^O Universe: Due to the recent exciting observational indications for 
a positive cosmological constant and the important consequences that this has 
to our understanding of the Cosmos, we will briefly present this model. 

Originally, the cosmological yl-parameter was introduced ad hoc by Einstein 
in his field equations in order to get a static solution (R = R = 0). From (||) 
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he derived R = (k/Ac) 1 / 2 and inserting this into ([?]) he obtained, in a matter 
dominated (p = 0) Universe: 

A c 2 

'-ra (24 > 

where /l c is the critical value of A for which R = R = 0. However, it was found 
that his solution was unstable and that small perturbations of the value of A c 
would change drastically the behaviour of R. From (0) we see that if k < 0, then 
R 2 is always nonnegative for A > 0, and thus the universe expands for ever, 
while if A < then the universe can expand and then recontract again (as in 
the k = 1, A = case). 

The recent SNIa observations (see section 3.2) and the CMB power-spectrum 
results (see section 3.1) have shown that the Standard Cosmological paradigm 
should be considered that of a flat, J? a — 0.7, f2 m ~ 0.3 model. Thus we will 
consider such a model in the following discussion. Evaluating (0) at the present 
epoch, changing successively variables: x = R 3 ^ 2 , y — x(Q m / ft^) 1 / 2 R 3 ^ 2 and 
6 = sinh -1 y and then integrating, we obtain: 



: sinh 



(25) 



and 



R — R 



sinh 3 



(26) 



It is interesting to note that in this model there is an epoch which corresponds 
to a value of R = R\ 1 where the expansion slows down and remains in a quasi- 
stationary phase for some time, expanding with R > thereafter (see Fig.l). At 
the quasi-stationary epoch, called the inflection point, we have R = and thus 
from (Q) by differentiation we have: 



Ri = 



2 Q h 



Rc 



(27) 



Now from (E5J) and (|27]) we have that the age of the universe at the inflection 
point is: 



3H y/f7X 



■. sinh 




(28) 



The Hubble function at t\ is: 



H{h) = iW^ACoth ( 3Ho ^ f, 



^ =>■ Hi — H \J 3i7^ 



so if t > t\ we must have H D < H\. 

This is an important result because it indicates that introducing an .!?A-term, 
and if we live at a time that fulfils the condition t a > tj, we can increase the age of 
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0.5 1 1.5 2 

t/t 

/ o 

Fig. 1. The expansion of the Universe in an Einstein de-Sitter (EdS) and in the pre- 
ferred A model. We indicate the inflection point beyond which the expansion acceler- 
ates. It is evident that in this model we live in the accelerated regime and thus the age 
of the Universe is larger than the Hubble time (H^ 1 ). 



the universe to comfortably fit the globular cluster ages while keeping the value 
of J2 m < 1 and also a flat (i?k = 0) space geometry. From (^8|) and (|23| ) and for 
the preferred values (1\ — 0.7 and Q m = 0.3 we indeed obtain t /t\ ~ 1.84 (see 
also Fig.l), which implies that we live in the accelerated phase of the Universe. 
Note that in order for the present time (t ) to be in the accelerated phase of the 
expansion we must have: fi\ > 1/3. 



Importance of k and A terms in global dynamics: Due to the recent 
interest in the A > 0, k ~ Universes, it is important to investigate the dynam- 
ical effects that this term may have in the evolution of the Universe and thus 
also in the structure formation processes (see Fig. 2). We realize these effects by 
inspecting the magnitudes of the two terms in the right hand side of (Q). We 
have the density term: 



SirGp 8nGpo 



and from (Effl) we have 



3 



(1 + zf = H^[2 m (1 + zf 



H 2 (l-f2 m ) (29) 



By equating the above two terms we can find the redshift at which they have 
equal contributions to the dynamics of the Universe. Evidently this happens only 
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Fig. 2. The strength of the three factors shaping the recent dynamics of the Universe. 
Compare the strength of the p and A term (k = 0) and of the p and k term. We have 
assumed H a = 72 km s _1 Mpc -1 , J? m = 0.3 and Qa = 0.7. 



in the very recent past: 

-1 (30) 



On 

Observations suggest that i? m ~ 0.3 and J?a — 0.7, and therefore we have 
z c ~ 0.3, which implies that the present dynamics of the universe are dominated 
by the yl-term, although for the largest part of the history of the Universe the 
determining factor in shaping its dynamical evolution is the matter content. 
Similar results are found for the fc-term in A = models. In this case we have 
from (H) that 

^ = i? 2 (l + 2) 2 |l-^ m | 

and thus the redshift at which the density and curvature terms have equal impact 
in the global dynamics, is: 

Zc = JT 2 - 

We see that as z increases the density term grows faster than the curvature term 
which is important only very recently. A similar line of argument shows that also 
in the radiation dominated era the A and k terms do not affect the dynamics of 
the Universe. 



Density parameter as a function of z: From ( |18[ ) and eliminating the curva- 
ture term, using (2\. = 1 — Q m — J?Aj we obtain the time evolution of the density 
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parameter J7 m : 

h 2 = Hi [rt m (i + zf + (i - « m - /2 A ) (i + zf + n A ] (3i) 

and using (|l5|) and ([!(]) we have: 

fl W = fl m + fl A /(l + z) 3 

i; r2 m + (i-A n )(i + z)- 1 + /2 A (i + z)- 1 [i/(i + 2) 2 -i] 1 j 

It is easy to see that whatever the value of Q m , at large z we always have 
Q{z) = 1. 

1.4 Distribution of Matter in the Universe 

The task of determining and quantifying the 3-dimensional distribution of extra- 
galactic objects, is not a trivial. Limitations are imposed on us by our position 
in space-time, and also the fact that we are located in the plane of a dusty spi- 
ral galaxy, which means that the light-absorbing interstellar dust will cause the 
magnitudes and sizes of galaxies to decrease, especially if we observe towards 
the Galactic plane. 

Generally in order to be able to quantify the geometry and topology of the 
large-scale structure of the Universe and to discriminate between the competing 
scenarios of structure formation, two at least issues should be addressed: 

(a) The quantification, in an objective manner, of the observed structure on large 
scales. For this we need to observe the positions and redshifts of galaxies tracing 
large- volumes, large enough to be considered a fair sample of the Universe. If 
not we can be influenced by local anomalies (like a local attractor or a local 
underdense region) and thus the interpretation of our results could lead us to 
erroneous conclusions. Today there are large samples available, complete to some 
apparent magnitude limit which are deep enough for volume-limited samples to 
be extracted from them (for example the PSCZ, ORS, LAS CAMPANAS, APM, SSRS, 
2dF, sdss catalogues). 

(b) The application of well defined and easily applicable statistics, which are able 
to distinguish between the different theoretical scenarios of structure formation 
and to compare between these scenarios and the observed topology and geometry 
of the large-scale structure of the Universe. 

Of historical interest are the first attempts to quantify the distribution of 
galaxies by Bok and Mowbray (in the mid 1930s). They compared the variance 
of galaxy counts N with that expected from a Poisson distribution (ie, compared 
(TV) 1 / 2 with o-=[jj J2( N i~ (^0) 2 ] 1/2 > where M is the number of fields) and con- 
cluded that a ^> {N) 1 / 2 an indication that galaxies cluster. More sophisticated 
statistics appear in the mid-50's when the large Lick galaxy counting project 



started giving its first results (cf. [114]). Through the years a wide variety of 



statistical techniques have been developed in order to study the geometry and 
topology of the distribution of matter on large-scales (for reviews see p2| [p3 ) . 
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Galaxies: Galaxies are the basic units of the distribution of matter in the 
universe. They are conglomerations of stars, dust and gas having typical sizes of 
a few kpc. Only recently, in the 1920s, was it realized that the fuzzy nebulae, that 
we now know to be galaxies, were not part of the Milky- Way. There appear to 
be many different types of galaxies, having different shapes, amounts of gas and 
dust, as well as star-forming regions. Many attempts have been made to produce 
a reliable classification scheme that bears relation to physical quantities. Galaxies 
come into three main categories: ellipticals, spirals and irregulars, in which the 
basic difference, as their names imply, is their apparent shape. However, further 
subdivision is essential in order to classify reliably the large variety of galaxies. 

The first attempt to classify galaxies was that of Hubble in 1926. His princi- 
pal morphological classification was provided by the "tuning-fork" diagram, in 
which the position of a galaxy was determined by the size of its nucleus and on 
the spiral arm tilt. For historical reasons, galaxies at the left of the diagram are 
called "early types" while at the right "late types" . Although this diagram was 
produced on the basis of appearance there seem to be a link between morphology 
and galaxy evolution since the integrated colours and spectral types monoton- 



ically increase from left to right ([178 and references therein). The fraction of 
the different galaxy types is: ~ 10% ellipticals (E), ~ 20% SOs, ~ 65% spirals 
(S) and ~ 5% irregulars (Irr). 

The spiral galaxies consist of three components: disk, bulge and halo. The 
disk is very thin (a few hundred pc thick). The surface brightness of stars in a 
typical disk falls exponentially with radius: 

I a = Jo cxp [-r/b] (33) 

where b is the scale-length of the disk (typically ~ 4 kpc). The spiral structure 
is thought to be the result of rotating density waves which produce shocks in 
the gas-rich disk which leads to star formation. 

The elliptical galaxies however, are ellipsoidal systems with an old population 
of stars and very little gas. They have a de Vaucouleurs luminosity distribution: 

I s = J e exp [-7.67((r/r ) a25 - 1)] (34) 

where r c is the effective radius enclosing half of the total light. Therefore J(r) 
falls off more slowly than r~ 2 for r < r and more rapidly for r > r c (this formula 
also fits the bulges of SOs and spirals). The mass of elliptical galaxies can vary 
widely; from dwarfs with M ~ 10 7 M Q to supergiants with M ~ 10 12 M 



Groups & Clusters of Galaxies: Galaxy clusters occupy a special position 
in the hierarchy of cosmic structures in many respects. Being the largest bound 
structures in the universe, they contain hundreds of galaxies, large amounts 
of Dark Matter (DM) and hot X-ray emitting gas, and thus can be detected 
at large redshifts. Therefore, they appear to be ideal tools for studying large- 
scale structure, testing theories of structure formation and extracting invaluable 
cosmological information (see recent reviews [24| , ||) Groups of galaxi es a re 
typically systems that contain < 10 — 20 galaxies (for a recent review see |105| |). 
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There are several classification schemes for clusters. A simple and clear-cut 
one is between regular and irregular clusters. 

® Regulars: They have smooth and symmetric structure, with high central 
density (> 10 3 /Mpc 3 ), a small fraction of spiral galaxies (< 20%), high veloc- 
ity dispersion (~ 10 3 km/sec) and high X-ray luminosity. The high velocity 
dispersion as well as the smooth structure is considered to be evidence that 
these clusters are in virial equilibrium (relaxed). About 50% of all clusters 
are regulars. 

(g) Irregulars: They have lumpy structure with evident substructures, lower 
velocity dispersion and lower X-ray luminosity. The fraction of spirals is 
higher than that in regular clusters (> 40%). 

Another distinct class of clusters is those having a central very bright galaxy 
(BCG or cD's). The BCG galaxies are giant ellipticals, some with multiple nu- 
clei (cD's) and masses M > 10 12 M©. There are different views regarding their 
formation history; one advocates that they form by "galactic cannibalism" where 
dynamical friction causes cluster galaxies to spiral towards the cluster core and 



then strong tidal effects take place (cf. [ 101 , U ) , the other advocates that they 



form at special locations, which are the kernels around which clusters will even- 



tually form by anisotropic accretion of matter (cf. [184 and references therein). 
Furthermore: 

♦ There is a well defined relationship between galaxy density and fraction of 
Hubble types. The fraction of E's and SO's increases smoothly with increas- 
ing galaxy-density ([Q). An interesting question is whether the above rela- 
tion is a direct result of the galaxy formation processes, or an evolutionary 
effect. 

♦ Clusters are highly flattened systems even more than elliptical galaxies. This 



fact is not due to rotation (cf. 149], [ 141 1 ) , and therefore it either reflects the 



initial conditions of their formation and/or the tidal effects acting between 
proto-clusters. In fact, the high peaks of an initial Gaussian density field, the 
possible sites of cluster formation, are non-spherical (J|]). The cluster mean 
projected ellipticity is e ~ 0.5 and there is evidence that clusters are more 
prolate-like (cf. @). 

Neighbouring clusters tend to point to each other up to distances of ~ 20—30 
h,- 1 Mpc (cf. Q). Dynamically young (irregular) clusters show a tendency 
to be more aligned with their neighbours and are preferentially found in 



high-density environments (]125|, [157]). 



For large-scale clustering, dynamical and cosmographical studies, it is ex- 
tremely important to compile large, whole-sky catalogues of groups and clusters. 
One of the first and extensively studied such catalogue is the ABELL / ACO cat- 
alogue Q, which was based on an eyeball selection procedure from sky survey 
plates. This all-sky sample contains 4073 rich clusters, nominally complete to a 
redshift z — 0.2. Its obvious limitations, due to the eyeball selection procedure, 
were superseded by the objectively selected APM clusters ([0) which were based 
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on the APM galaxy catalogue containing 4 x 10 6 galaxies (p^j). This cluster cata- 
logue covers latitudes b < —35° and contains 950 clusters, typically poorer than 
the ABELl/aco one's. Furthermore, with the Rosat whole-sky X-ray survey it 
was possible to construct X-ray selected cluster catalogues (cf. Q), which suffer 
less from projection effects, which can produce phantom clusters. For example 
the Reflex X-ray cluster sample contains ~ 450 clusters (pl|). 



Fig. 3. (Included separately as a JPG image) An example of a filamentary supercluster: 
The 2D distribution of APM clusters members. The inserted plots show the smooth 
2D galaxy density map of some of the cluster members. For the case of A3128 we 
overlay the X-ray emission contour plots. The cluster's elongation along the filamentary 
supercluster is evident. 



Superclusters, Filaments & Voids: Superclusters are aggregates of clusters, 
groups and galaxies. They are the largest, isolated, but dynamically un-relaxed 
due to their size, objects in the large scale distribution of matter and thus they 
are ideal probes of the initial conditions that gave rise to them. This fact is 
because typical peculiar velocities of clusters are v pec ~ 10 3 km/sec and therefore 
in a Hubble time (1/H a ) they can move no more than ~10 hr 1 Mpc, which is 
substantially smaller than the scale of a typical supercluster. Superclusters can 
be identified in 3-D catalogues of clusters but also in 2-D projections of galaxy 
distributions. Regions of high density in clusters can be identified also on the 
APM galaxy map. 

The large scale clustering pattern of galaxies is expected to be character- 



ized by a filamentary and sheet-like distribution (cf. [194 , |p5|, |17£|). Indeed 



many authors have been finding that the vast majority of the superclusters are 
flattened with a predominance of filament- like shapes (cf. 13 and references 
therein). Figure 3 shows the 2-dimensional projection of a filamentary super- 
cluster containing 5 APM clusters. 

Superclusters are not centrally condensed objects (like clusters) and their 
typical size is ~ 30 — 50 h^ 1 M pc. H owever larger structures, with a length 
~ 200 h^ 1 Mpc, may exist (cf. ||l76|] and references therein). Detailed stud- 
ies have shown that elongated bridges (~ 30 hr 1 Mpc) of galaxies have been 
found to connect rich clusters. Since wide-angle three-dimensional surveys be- 
came available, the filamentary distribution of galaxies has been a constantly 
observed feature. Even the original Cfa survey (cf. j72|) showed netwo rks of 
filaments mostly connecting rich clusters of galaxies and large voids (cf . |116|| ) , 
a fact which has been confirmed by all the recent surveys (SSRS, ORS, PSCZ etc). 
Voids are regions of density well below the average value. In all deep radial- 
velocity surveys, the velocity distribution shows striking empty regions were no 
(or very few) galaxies exist. It is an extremely difficult task to identify voids in 
2-D projections of galaxy distributions since the projection will tend to smooth- 
out any such structure. An extremely interesting question, relevant to theories 
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of structure formation, is whether the voids are empty of luminous matter or 
empty of all matter. 
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2 Distance Scale, Hubble Constant & the Age of the 
Universe 

One of the most important parameters in determining the fate of the Universe as 
a whole, is the present day expansion rate of the Universe, which is encapsulated 
in the value of the Hubble constant H . Its value sets the age of the Universe 
and specifies the value of the critical density p CT , and through this route the 
geometry of the Universe. 

From the Hubble law (||) it is evident that in order to determine its value 
we need to determine the expansion velocity (redshift) as well as the distance 
of extragalactic objects, but within a redshift such that space-curvature effects 
do not affect distances (^3|). A further concern is that local gravitational effects 
produce peculiar velocities, that are superimposed on the general expansion. This 
can easily be seen in the toy model of (|2]) , if we allow \ = x{t)- Then 



with the factor on the right being the peculiar velocity. Since the observer as 
well as the extragalactic object, of which the distance we want to measure, have 
peculiar velocities then the above equation becomes: 



where v(0) is the velocity of the observer and f is the unit vector along the line- 
of-sight to the extragalactic object. It is then obvious that in order to measure 
H , the local velocity field should be measured and the extragalactic distances 
corrected accordingly. It is easily seen that if both observer and galaxy take 
part in a coherent bulk flow, having the same amplitude at the observer and 
galaxy positions, then the right-hand part of ( |35| ) vanishes. In general however, 
one needs good knowledge of the velocity field in order to correct distances 
adequately. 

2.1 Distances of Extragalactic Objects 

Our only means of obtaining information and therefore knowledge of the struc- 
ture and dynamics of the Universe (on all different scales) are through the elec- 
tromagnetic radiation that we receive. Therefore it is of primary importance 
to define a system of measuring luminosities taking also into account that the 
Universe expands and that light loses energy through a variety of processes. 

If we assume that light propagates with no loss of energy, then the apparent 
luminosity of a source I, is related to its absolute luminosity L, by: 



v = H d + Rx , 



cz = H o d+(v(d)-v(0))-r 



(35) 



I = 



L 



(36) 



47rr 2 



where r is the distance to the source. We can see the extreme importance of 
determining the pair (/, L), since such knowledge would provide the distance of 
the source, r. 
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Due to historical mostly reasons we use a logarithmic brightness system by 
which an object with a magnitude of 1 is 100 times brighter than an object with 
a magnitude 6. We have: 



m = —2.5 log 10 I + ci 



M 



2.5 log 10 L + c 2 



(37) 



where to is the apparent magnitude and M the absolute one. Therefore, using 
( |36| ) we have: 

m - M = 5 log 10 r + c 3 (38) 

where c\ t % are constants which depend on the filter used and C3 is that value 
for which m — M at a distance of 10 parsecs (see section 2.3) from the Earth, 
and thus C3 = —5. In extragalactic astronomy, instead of pc, we use Mpc and 



therefore (38) becomes: 



M = 5 log 10 r + 25 



(39) 



The above definitions are somewhat ideal, since in the real world we do not 
observe the total apparent magnitude, but that corresponding to the particular 
range of spectral frequencies, that our detector is sensitive to, and those allowed 
to pass through Earth's atmosphere. If we detect l v over a range of frequencies 
(v ± Sv), then the observed apparent magnitude is m — —2.5 log 10 J g l v Av + c. 
However, neither the atmosphere nor the detectors have a sharp v limit and 
therefore it is better to model these effects by a sensitivity mask F v , and the 
observed apparent magnitude is then: 



TO F„ 



-2.5 log 



10 



F v l v &v 



+ c 



(40) 



If F v = 1 then the apparent (or absolute) magnitude of a source is called the 
bolometric magnitude. 

How do the above definitions change by taking into account the fact that 
the Universe expands? To answer this we need a metric of space-time, which in 
our case is the Robertson- Walker metric. Since light travels along null-geodesies, 
a fundamental concept of distance can be defined by the corresponding light- 
travel time, which is called proper distance. If a light signal is emitted at a 
galaxy Q\ from the coordinate position (ri,8 0} 4>o) at time t — and received by 
an observer at Go at (ro, 6q,4>o)j then these events are connected only by the light 
signal and since all observers must measure the same speed of light, it defines a 
very fundamental concept of distance. Obviously, it depends on the curvature of 
space, and since ds = we have from (Fy): 



r'i 



( R(t) sin 



[ R{t) sinh~ ] 




(41) 



In the expanding Universe framework, the expressions (36) and (|3§|) change, to: 

L 



An dl, a (l + zY ' 



to — M = 5 log 10 g? l + 25 



(42) 
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where g?l = d plo (l 4- z) is the luminosity distance. It is obvious that the dis- 
tance measure of an extragalactic object depends on the underlying Cosmology. 



A proper derivation of the luminosity distance (cf. [ 183 1 , [ 1 12 1 ) provides the 
following expression: 



2c 



H Q 2 



Qz + {n- 2) [Vl + zQ - 1 



(43) 



where Q contains all the contributions (mass, energy density, curvature). 

Another important distance definition is that of the angular-diameter. It is 
based on the fact that a length I, subtends a smaller angle, 9, the further away 
it is (I oc 1/9) and it is given from: relation to: 



do 



(i + z y 



(44) 



Note that this notion of distance is used to derive the CMB power spectrum 
predictions of the different cosmological models (see section 3.1). 



2.2 Biases affecting Distance Determinations 

Many effects, related either to our position in space (being on the plane of a 
dust-rich spiral galaxy), to natural limitations (introduced for example by the 
expansion of the Universe) or to detector related issues, introduce systematic 
biases that affect our ability to measure accurately distances. Below, I list a few 
of the most important such effects. 



K-correction: Since bolometric magnitudes are not possible to measure, but 
rather magnitudes over a particular wavelength range, it is important to correct 
these magnitudes for the effect of the expansion of the Universe. These con- 
siderations result in modifying the distance modulus by a factor, the so-called 
A'— correction factor: 

to - M = 5 log 10 d L + 25 + K{z) (45) 

This factor arises from the fact that when we measure the magnitude of sources 
at large distances and at a particular frequency, say v , we receive light emitted 
from a different part of the spectrum, v c . It could well be that in this latter part 
of the spectrum the extragalactic object is particularly fainter or brighter than in 
the nominal one, v a . Furthermore, a combination of different factors; evolution, 
intervening absorption processes or detector sensitivity for example, result in 
energy losses as a function of wavelength, which can be expressed by a mask 
F(v ) (fief). Knowing F(v ) one can estimate the A-factor by integrating the 
spectrum at the source rest frame. For example, such calculations have shown 
that a typical value for spiral galaxies at z = 1 is A 2 (in general K(z) cx z). 
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Malmquist bias: Due to the nature of astronomical observation there is a 
minimum flux density above which we select extragalactic objects. As we look 
deeper in the Universe we tend to pick up a relatively larger fraction of intrin- 
sically bright objects (ie., only the brighter end of the luminosity function is 
sampled). This bias arises when determining distances from apparent magni- 
tude limited samples. If the individual absolute magnitudes Mi of a sample of 
extragalactic objects have a Gaussian distribution around (M) with dispersion 
a, then this bias, for the case where the distribution of extragalactic objects is 
considered homogeneous, is given by: 

A(M) = 1.382er 2 (46) 

How does this bias affect the determination of extragalactic distances? The in- 
ferred distances of extragalactic objects are typically smaller than their true 
distances. From (|3^ ) we have that: 

r ~ r inl.382<T 2 /5 
' cor ~ ' raw ±lJ 

We illustrate this bias in Fig. 4; as the distance increases, M(rnu m ) becomes 
brighter and therefore the brighter end of the luminosity function is sampled. 
For a larger mn m (deeper sample) the value of M(mn m ) increases (less luminous) 
so we have a smaller A(M). Conversely, for a given m\ lm and (M), the bias 
increases with distance. 
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Note that we have considered a fairly straight-forward case, ie. that of a 
sample with a unique (M) value. In real samples of extragalactic objects we 
have a range of such values and therefore this bias is not easily seen. 

A related bias that also affects extragalactic distance determinations is the 
fact that there are larger numbers of objects at larger distances and therefore 
within a given range of estimated distances, more will be scattered by errors 
from larger to smaller distances than from smaller to larger ones. 

Galactic absorption: Interstellar gas and dust absorbs the background light 
with dust scattering more efficiently the blue light, and thus the background 
light appears artificially reddened. From simple geometrical considerations it is 
easy to show that the flux l v of an extragalactic source, transversing a Galactic 
layer of thickness ds, at an angle b from the equatorial plane, suffers losses 
Sl v /l v oc ds cosec(fr) and therefore: 

^p- — — /t„£„cosec(&) (47) 

where the constant of proportionality k v is the absorption coefficient at the 
spectral frequency v. Therefore, integrating we have: 

l v = 1° exp [-A cosec|6|] (48) 

where the integration constant 1° is the incident and l v is the observed flux, 
while A — J n v ds is the optical thickness. Therefore to take into account this 
effect ( p6| ) should change to: 

h - = exp [-Icosec|fr|] = CXP [A C ° SeC|6|] • (49) 

Values of A slightly vary for different spectral frequency bands, but a generally 
accepted value, in V, is ~ 0.2. We see from ( fil] ) and ( |39| ) that 

^*truc ~ ^raw exp [— A cosec|6|/2] (50) 

ie., the distance of an extragalactic source at a given absolute magnitude can 
be significantly overestimated at low galactic latitudes if this effect is not taken 
into account. Note however that the cosec(fr) model is oversimplified since the 
distribution of gas and dust in the Galaxy is rather patchy. 

Cosmological Evolution: As we look back in time we see a distribution of 
extragalactic objects (normal galaxies, AGN's, clusters) in different evolutionary 
stages. It may well be that their luminosity and/or mean number density is a 
function of cosmic time, a fact that will affect distance determinations based on 
local calibrations of the relevant scaling-relations. 
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Aperture effect: Since galaxies are extended objects with no sharp outer 
boundaries, their photometric measures will depend also on the size of the tele- 
scope aperture since at different distances different fraction of a galaxy will fit 
in the aperture. This is a distance-dependent effect, since diameters scale like 
1/r, and therefore it may affect the distance estimate. 



2.3 Distance Indicators: 



In order to develop the distance ladder from local to distant cosmic objects, one 
starts from the local distance scale (for a detailed account see [143] [144|). 



Galactic Distance Indicators: The primary method used to estimate the 
distances to nearby stars is that of Trigonometric Parallax. As the Earth orbits 
the Sun, we view the Universe from different points of the orbit throughout 
the year. This difference is maximum every 6 months when the Earth is on 
the opposite side of its orbit around the Sun. This kinematic method provides 
the basic unit of distance, the parsec, used by astronomers. It is defined as the 
distance at which a star would have a trigonometric parallax of 1 (parsec = 
parallax + second): 

lpc = 3.086 x 10 13 km = 3.26 light years 

Note that this method is effective out to ~ 60 parsecs and that the nearest star 
to us (ce-Centauri) is at a distance of 0.75 parsecs. 

Among the many distance indicators, used to determine distances within 
our Galaxy, a particularly interesting one is the Main sequence fitting method. 
This takes advantage of the fact that stars in globular clusters are at a common 
distance and that there is a unique correlation between spectral stellar type and 
absolute luminosity (the H-R diagram). Therefore by measuring the distance, via 
a kinematic method, to one relatively near globular cluster, one sets the zero- 
point of this method and then by observing the apparent magnitude - spectral 
type distribution of other globular clusters, one can determine their distance. 



Extragalactic Distance Indicators: The next step is based on Cepheid Vari- 
able Stars. This method has been traditionally used within our Galaxy and in 
the nearby Large Magellanic Cloud (LMC), but with the Hubble Space Telescope 
it has been successfully used out to ~ 20 Mpc (cf. @|). A strong and tight re- 
lationship exists between the intrinsic luminosity of these stars and their period 
of variation (pulsation) which results in a Period-Luminosity relation: 

L cx log P 13 =>■ log L = zero-point + 1.3 log P (51) 

Once this relation has been calibrated, it provides the absolute luminosity of the 
distant Cepheid stars and via (p9|) , the distance of their host galaxy. Although 
this relation has a scatter, in the /—band, of only ±0.1 mag, systematic effects 
may exist. For example, a serious concern is whether there is any environmental 
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dependence of the relation. It has been suggested that a different metalicity of 
the host galaxy may significantly affect the zero-point of the relation and thus the 
determined distance. These effects can be taken into account and this method has 
proved to be fundamental in the recent determinations of H , because it provides 
the link between the primary galactic indicators and the local extragalactic ones, 
which then provide the calibration for other secondary indicators operating in 



much larger distances (cf. |61| and references therein). 

In developing the distance scale, we now need effective indicators that can be 
used to very large distances. Other scaling relations have been found between a 
distance dependent (ex. brightness, diameter) and a distance independent (ex. 
rotational velocity, stellar velocity dispersion) quantity. It is evident that from 
such relations one can extract distance information. The main assumption in 
such a use of these scaling relations is that they are not environment-dependent 
(which has been shown to be a valid assumption). Such a relation for spiral 



galaxies is the Tully- Fisher relation [175|, which relates the rotational velocity 



of a spiral to its total infrared luminosity: L„ cx V^ ot or its total blue luminosity: 
Lb oc V T " t with a ~ 2.4 — 2.8. It has a reasonable theoretical justification: 
Rotational velocities in spirals are related to mass according to 

K 2 ot cc |- (52) 
Assuming that all spirals have the same surface brightness S, then 

S cx L/R 2 LozR 2 (53) 
If the mass to light ratio is constant then 

^cxK 2 ot£ 1/2 ^«K 4 ot (54) 

For ellipticals a similar relation holds, the Faber-, Jackson relation p(i| |, which 
relates the absolute luminosity of the galaxy with the stellar velocity dispersion 
(L cx cr 3 ~ 4 ) or a variant, the so-called D n — a relation (cf. pofl), which relates 
the diameter D n of an elliptical (defined as that within which the mean surface 
brightness is 20.75 mag arc sec~ 2 in B) to the stellar velocity dispersion a: 
D n oc a x with x ~ 1.2—1.3 . The typical accuracy of these distance estimations 
is ~ 20% and the usual assumption is that they do not evolve with redshift over 
the scales used, and that systematic effects can be corrected for. However, there 
are some indications for an evolution of the -B-band Tully-Fisher relation in the 



interval studied 0.1 < z < 1 |195] 



Another very important distance indicator, which can be used to large dis- 
tances, is the Supernova SNIa brightness at maximum. It is thought that the 
maximum luminosity of such supernovae (explosion of white dwarfs in binary 
systems which become gravitationally unstable - reaching the Chandrasekhar 
limit - due to the accretion of matter from the secondary - see j9^] ) is a Univer- 
sal constant and since the intrinsic luminosity of a SNIa is high, they can be seen 
out to cosmological distances. Furthermore, a correlation was found between the 
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Fig. 5. The SN la Hubble diagram vs Cosmological model predictions (from [HE] with 
permission) . 



supernova luminosity and the brightness decay time, which provides a further 
parameter that reduces the scatter in luminosity to ±0.3 mags. Using this dis- 
tance indicator one can construct the Hubble diagram to very large-distances 
and thus determine the deceleration of the Universe by mapping the region of 
the Hubble diagram that deviates from linearity. This has been recently achieved 
by two independent groups, the Supernova Cosmology Project - SCP [HE] and 
the High-z Supernova Search Team - HZT [156], which have found that the de- 
rived Hubble diagram is that expected from a accelerating expansion, which can 
be provided by a non-zero cosmological constant (see Fig. 5). 

Although I do not plan to present all the secondary distance indicators, one 
that is potentially very important and susceptible to small uncertainties and 
systematics, is the surface brightness fluctuation method. This method is based 
on the fact that the discreteness of stars within galaxies depends on distance. 
This methods has an accuracy of ~ 5% in distance. 



Direct Distance Indicators: Clusters of galaxies are filled with hot and trans- 
parent gas which can be fitted by a thermal Bremsstrahlung spectrum with 
T~ 5 x 10 7 — 10 8 K (covering a range of 2 - 8 KeV). The physics of the hot 
intcrcluster gas provides the means of measuring directly the distance of clus- 
ters, without need of intermediate steps. This method is based on the so called 
Sunyaev-Zeldovich effect which is the distortion of the original CMB spectrum, 
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by the Compton scattering to higher energies of the CMB photons from the 



hot electrons of the plasma (for recent reviews see [ 13 1 1 , |L7]]). This distortion 
decreases the brightness of the CMB spectrum at the longer wavelength range 
while it increases the photon energies in the shorter wavelength range. Lets as- 
sume a cluster of radius R at a distance D, subtending an angular separation 9 
on the sky. The change of the CMB brightness temperature is proportional to 
the line integral of electron number density through the cluster: 

AT f 

— oc - J T e {l) n e d (55) 

then from isothermality we get AT jT oc T c n c R. Furthermore the luminosity 
of the Bremsstrahlung radiation together with the assumption of isothermality, 
gives: 

From (pq), the observed X-ray flux (F K oc L x /D 2 ) and R = 9D we have: 

F^ntT^8*D*^ 2 -±- (57) 



and solving for D we get: 

D oc 



I AT 

{ — J T 3/2 



Therefore measuring F x , 6, T c we obtain an absolute determination of the cluster 
distance. The temperature, T e , can be measured either from the shape of the X- 
ray continuum or from line emission (especially of iron). Note however, that we 
have assumed sphericity, isothermality and a smooth distribution of n c . Most 
clusters are flattened (cf. Q), and show significant substructure apparent in 
the optical and X-ray images. In such clusters the above procedure may provide 
highly uncertain and biased distance estimates (cf. J77J). A recent study of ~ 100 
clusters has shown that once cooling flows are taken into account, isothermal 



profiles fit well ~ 90% of the clusters ([185 , but see also |4 



2.4 The Value of H Q and the Age of the Universe 

The immense effort that has been put towards the goal of determining H , with 
an accuracy of a few %, has recently come to fruition with different methods 
providing consistent (within < 10%) values. In Table 1, I present a list of some 
of the most recent determinations of H . Note that there are 3 different mea- 
surements based on SNIa, giving different values of H . Although many of the 
SNIa they use are common, the difference is most probably attributable to the 
different local calibrations that they employ. Thus, the differences in their de- 
rived -ffo's should reflect the systematic uncertainty introduced by the different 
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local calibrations and it is indeed comparable to the systematic uncertainty that 
the individual studies have estimated. 

The method in the last row of Table 1 is based on gravitational lensing 
(cf. ^0|). The basic principle behind this method is that there is a difference 
in the light travel time along two distinct rays from a source, which has been 
gravitationally lensed by some intervening mass. The relative time delay (.4t), 
between two images of the source, can be measured if the source is variable. 
Then it can be shown that the Hubble constant is just: 

Aff 1 

H =C— (58) 

where A6 is the image separation and C is a constant that depends on the lens 
model. Although, this method has well understood physical principles, still the 
details of the lensing model provide quite large uncertainties in the derived H Q . 

A crude iV-weighted average of the different iJo-determinations in Table 1 
gives: 

H = 68 ± 6 km sec" 1 Mpc" 1 

where the uncertainty reflects that of the weighted mean (the individual uncer- 
tainties have not been taken into account). However, there seems to be some 
clustering around two preferred values (H ~ 60 and ~ 72 km s" 1 Mpc -1 ) and 
thus the above averaging provides biased results. More appropriate is to quote 
the median value and the 95% confidence limits: 

H = 72^ km sec" 1 Mpc" 1 , (59) 

The anisotropic errors reflect the non-Gaussian nature of the distribution of the 
derived if -values. Note that the largest part of the individual uncertainties, 
presented in Table 1, of all methods except the last two, are systematic because 
they rely on local calibrators (like the distance to the LMC), which then implies 
that a systematic offset of the local zero-point will "perpetuate" to the secondary 
indicators although internally they may be self-consistent. A further source of 
systematic errors is the peculiar velocity model, used to correct the derived 



distances, which can easily introduce ~ 7% shifts in the derived H values [146 
With the value (p9), we obtain a Hubble time, £h, equal to: 



ta = H- 1 = 13.6±g. 6 Gyr's (60) 

the uncertainties reflecting the 95% confidence interval. 

It is trivial to state that the present age of the Universe t , should always 
be larger than the age of any extragalactic object. A well known problem that 
has troubled cosmologists, is the fact that the predicted age of the Universe, in 
the classical Einstein de-Sitter Cosmological Model is smaller than the measured 
age of the oldest globular clusters of the Galaxy. This can be clearly seen from 
© and 

e* S = k^. 7 Gyr's (61) 
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Table 1. Some recent determinations of the Hubble constant, based on different meth- 
ods. 



Method 


N 


H 


(*) 


reference 


Cepheid 


23 


75±10 


0.006 


Freedman et al 2001 


2-ary methods 


77 


72±8.0 


< o.i 


Freedman et al 2001 


IR SBF 


16 


76±6.2 


0.020 


Jensen et al. 2000 


SN la 


36 


73±7.3 


< o.i 


Gibson & Brook 2001 


SN la 


35 


59 ±6 


< o.i 


Parodi et al 2000 


SN la 


46 


64 ±7 


< o.i 


Jua et al 1999 


CO-line T-F 


36 


60±10 


< 0.11 


Tutui et al 2001 


S-Z 


7 


66±15 


< 0.1 


Mason et al 2001 


Grav. Lens 


5 


68 ± 13 




Koopmans & Fassnacht 2000 



and although the latest estimates of the globular cluster ages have been drasti- 
cally decreased to Q: 

t gc ~ 12.5 with 95% lower/upper limit is 11/16 Gyr's . 

One should then add the age of the formation of the globular clusters and as- 
suming a redshift of formation z ~ 5 then this age is ~ 0.6 — 0.8 Gyr's which 
brings the lower 95% limit of i gc to ~ 11.6 Gyr's (see however |jo) for possible 
formation at z > 10). It is evident that there is a discrepancy between t and 
t gc . This discrepancy could however be bridged if one is willing to push in the 
right direction the 95% limit of both % and t gc . 

However, other lines of research point towards the age-problem. For example, 
if at some large redshift we observe galaxies with old stellar populations, for 
which we know the necessary time for evolution to their locally "present" state, 
then we can deduce again the age of the Universe. In an EdS we have R cx t 2 / 3 
and thus we have: 

t = t z (l + z) 3/2 (62) 

Galaxies have been found at z ~ 3 with spectra that correspond to a stellar 
component as old as ~ 1.5 Gyr's, in their local rest-frame. From (|6^) we then 
have that t ~ 12 Gyr's, in disagreement with the EdS age (|6l|). 

This controversy could be solved in a number of ways, some of which are: 

• invoking an open (J? < 1) cosmological model, 

• assuming that we live in a local underdensity of an EdS Universe, 

• invoking a flat model with Q\ > 0. 

The first possibility is in contradiction with many observational data and 
most importantly with the recent CMB experiment s (boome rang, maxima 
and DASi), which show that fi = 1 (see @ 



0| |6j (13 



The second possibility [173| can solve the age-problem by assuming that we 
live in a local underdense region that extends to quite a large distance, which 
would then imply that the measured local Hubble constant is an overestimate 



28 M. Plionis 



of the global one by a factor: 

sh a°- 6 SN 



H 36 N 

where the bias factor 6, is the ratio of the fluctuations in galaxies and mass. To 
reduce the Hubble constant to a comfortable value to solve the age problem, say 
from 72 to 50 km sec -1 Mpc -1 , one then needs SN/N ~ —0.9/6. Values of 6 are 
highly uncertain and model dependent, but most recent studies point to 6 ~ 1 
(cf. [p8| ) which would then mean that we need to live in a local very underdense 
region, something that is not supported by the linearity of the Hubble relation 



out to z < 0.03 (cf. BJ) or out to z < 0.1 (cf. 1 171 ). This is not to say that we 



are not possibly located in an underdense region, but rather that this cannot be 



the sole cause of the age-problem [193 



Thus we are left with the last possibility of a Universe dominated by vacuum 
energy (a Universe with SI a > - see section 1.3). If we live in the accelerated 
phase (see Fig!) we will measure: 

H(t ) > H(t) with t > t 

ie., a larger Hubble constant, and thus smaller Hubble time as we progress in 
time, resolving the age-problem. In fact we have strong indications (see next 
section), from the SNIa results, which trace the Hubble relation at very large 
distances (see section 3.2), and from the combined analysis of CMB anisotropy 
and galaxy clustering measurements in the 2dF galaxy redshift survey (54) , for a 
flat Universe with J?a — 0.7. Then from ( |23| ) we obtain the age of the Universe 
in such a model: 

t£ = 1.446 x t™ s = 13.1±?;f Gyr's 

Indeed the resolution of the age-problem gives further support to the J?a > 
paradigm. 
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3 Determination of the Matter/Energy Density of the 
Universe 

The existence of large amounts of Dark Matter in the universe, manifesting 
itself through its gravitational effects, is a well established fact, although the 
precise amount has been a matter of a lively debate through the years. Attempts 
to identify the DM as normal baryonic matter has failed, mostly due to the 
nucleosynthesis constraints imposed by the successful hot Big-Bang model and 
the large temperature fluctuations of the CMB that it predicts in a flat Universe 
~ \^f)- Possibly some of the DM could be neutral hydrogen, in the form 
of Lyman-a clouds, but it is estimated that it could contribute only fl < 0.01. 
Similarly, the possible solid form of baryonic material (eg. dust grains, Jupiters, 
dwarfs with M < 0.08M Q or neutron stars) would contribute little to fl. 

Two recent determinations of the deuterium abundance, which combined 
with the BBN (Big-Bang Nucleosynthesis) predictions, give the total baryonic 



mass in the universe, have provided slightly discrepant results (see reviews [107], 
[ |163| ), covering the range: 

0.005 < f2 B h 2 < 0.024 (63) 

and therefore for h = 0.72 we have: 0.01 < J?b < 0.046. However, the recent 
analysis of the results from the BOOMERANG CMB experiment have provided a 
value mostly compatible with the lower deuterium abundance and thus higher 
fifth 2 value (see further below). 

In this section we present a variety of methods used to estimate either the 
total mass/energy density of the Universe or its mass density, i? m . 



3.1 The CMB fluctuation spectrum 

The most straight forward approach to estimate the total matter and energy 
density of the Universe (ie., the total J?) is by means of the measurement of the 
fluctuation spectrum of the CMB. Before recombination at z < 1100, the baryons 
and photons are tightly coupled, oscillating acoustically due to gravity (on sub- 
horizon scales). Only after recombination do the acoustic oscillations stop and 
the density fluctuations grow. The fluctuations emerging from the last scattering 



surface are a series of peaks and troughs [169 and as the different wave-lengths 
are projected to different angular scales on the last scattering surface and de- 
pending on the underlying cosmological model, they produce a characteristic 
structure of peaks on the CMB power spectrum (for a recent review see [Q 
and references therein). This method is in effect based in measuring the angular 
extent of a physical scale on the last scattering surface. The curvature of space 
enters through the angular distance to the last scattering surface. Therefore, the 
same physical scale will be projected to a smaller angular scale on the CMB sky 
in a positively curved background, while it will be projected to a larger angular 
scale in a flat or to an even larger scale in a negatively curved background space. 
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Fig. 6. CMB spectrum from the bommerang, maxima and dasi experiments with the 
error boxes of the measurements. The predictions of the popular inflationary model 
and one non-Gaussian (global texture) model (from b8[ with permission). 



To define the CMB power spectrum one starts by expanding the temperature 
fluctuations of the CMB sky in spherical harmonics: 

fry, oo m—£ 

— {6,4>)=Y J E OArJWM) > (64) 

1=1 m =-£ 

then if the fluctuations are Gaussian, the 2-point correlation function contains 
all the statistical information, and can be defined as: 

= ^(2£ + 1) Ci W e P t (8, 0) . (65) 

i 

where Wi is the window function representing the beam characteristics of the 
experimental apparatus used to observe the CMB sky, while the average is over 
all positions on the sky. One then invokes the ergotic theorem, ie, that the above 
average is equivalent to being over different realizations of our Universe. Then 
assuming random phases one can define the CMB power spectrum Ci as the 
ensemble average of the coefficients a,£ m : 

Q = (|a, m | 2 ) 
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Fig. 7. Different combinations of the cosmological 
CMB power-spectrum - degeneracy problem (form 



>arameters can result in the same 
99[ with permission) 



The different cosmological parameters will reflect onto a different structure 
of peaks in the structure of the CMB power spectrum. The position of the first 
peak is determined by the global mass/energy density of the Universe and the 
dependence of ^ pea k on Q can be approximated by: 



''peak 



220 
~7fl 



(66) 



Note however, that this approximation is not correct in yl-dominated universes 
and small corrections should be applied (cf. |99]). Many recent experiments like 
the bommerang, maxima and dasi (cf. Q, [43|, Jl64| (13§) find: 



fltot - 1-02 (±0.07) 

Many other cosmological parameters (for example f2 m , fl^ H a , baryon con- 
tent of the universe, the spectral index n of the inflationary perturbation spec- 
trum, etc) affect the structure of the peaks, beyond the first one (cf. j7lj). Deter- 
mining the CMB spectrum up to a few thousand i's can put strong constraints on 
these parameters. Current experiments trace the CMB spectrum up to I ~ 1000 
and indeed they have detected two more significant peaks at roughly I ~ 540 
and 840 E| (see Fig.6). 
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Note however, that different combinations of the cosmological parameters 
can conspire to produce exactly the same CMB spectrum; this is the so called 
degeneracy problem (see Fig. 7) and therefore in order to provide limits to these 
cosmological parameters one needs to assume priors and/or constrain different 
combinations of these parameters. However, the more accurate the derived CMB 
spectrum the weaker the necessary priors^. 

The latest data and CMB spectrum analysis provides very stringent con- 
straints to the baryon content of the Universe: n^ti 2 ~ 0.022^q'qq|, consistent 
with the primordial nucleosynthesis constraints (see |6^), and to the spectral 
index of the power spectrum of primordial perturbations: n ~ 0.96 ± 0.1 
Furthermore, combined analyses with other cosmological data, can be used to 
break the above mentioned degenerecies (see below). 



3.2 The Hubble diagram with SNIa 

As we have already discussed in section 2.3, the Hubble diagram of supernovae 
SNIa can be used not only to determine the Hubble constant (at relatively low 
rcdshifts) but also to trace the curvature of the Hubble relation at high redshifts 
(see |fL36| and references therein). 



2 With the new CMB experiments - MAP and PLANCK - the CMB power spectrum 
will be determined to an unprecedent d etail , providing extremely accurate values for 



more than 10 cosmological parameters [160] 
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The two groups working laboriously on this subject (SCP and HZT) have 
found consistent results, by which the distant SNIa's are dimmer on average by 
0.2 mag than what expected in a flat EdS model, which translates in them being 



10% further away than expected ([119|, |135| ). This implies that we live in 



an accelerating phase of the expansion of the Universe, a fact that supports a 
non-zero cosmological constant. The confidence intervals that their results put 
in the f2 m — Q\ plane are shown in Fig. 8. These results can be quantified by the 
following expression [|ll9| : 



0.8^ m - 0.6/2 A = 0.2 (±0.1) 

Together with the CMB fluctuation spectrum results we obtain: 

Q m ~ 0.29 H A ~ 0.71 . 

However, since our understanding of the physics of SNIa'a is not complete (cf. 
[ j92| , |pi| ) there could be some systematic effect, correlated with distance (eg. 
evolution), which could explain the dimming of the distant SNIa's and thus 
alleviate the J?a > interpretation. In Fig. 9 we show the distance modulus 
residuals after subtracting an open J? m = 0.3 Hubble relation. The systematic 
distant-dependent effect mimics the accelerated expansion Hubble relation out 
to z ~ 0.8 — 1. Beyond z ~ 1 the two relations depart due to the fact that 
the accelerated phase has to first pass from a decelerating one (see discussion in 
section 1.3) and this could provide a strong test for the possible distant depen- 
dent systematics. In fact, the recent discovery of the furthest known supernova 



(SN 1997ff) at a redshift of z ~ 1.7 [137|, has provided evidence of the decel- 



erating phase of the presently accelerating Universe (however, more very high-z 
supernovae are necessary to confirm this extraordinary result). 
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3.3 Clustering of Galaxies, Clusters &; QSO's 

If we assume a continuous density field of extragalactic objects, p(r), with mean 
p, we define the fluctuations of the field at position x as: 

*(*) = P -^I , ( 67 ) 
P 

Obviously we have that (S(x)) = 0. The correlation function is defined as: 

ar) = (6( X )5( X + r))= ^ X + r ^-^ . (68) 

P 

and quantifies the extend to which the density fluctuations at a given point are 
correlated to those at a distance r. The value of 2-point correlation function at 
zero-lag is therefore the variance of the random process: 

e(o) = (s 2 (x)) 

which measures the excursions of the density about its mean value. However, this 
is not a well defined quantity because usually the fluctuation field is smoothed 
to some resolution, say R. We then evaluate the variance of this field as r — ► 
(see fn] below). 

In many problems it is convenient to work in wave-number space. The Fourier 
transform of S(x) is: 

5 k = J 5{x)e tkx d 3 x (69) 
and it is convenient to separate S k in modulus and argument: 

S k = \S k \e^ , 

where e k are the phases (usually assumed to be randomly distributed in [0,7r), 
although the non-linear evolution of structure introduces phase correlations, cf. 
[|33|). The variance of the amplitudes is the power spectrum: 

P(k) = (|4| 2 ) 

which is the Fourier transform of the correlation function ( Wiener- Khinchin 
theorem) : 

P(k) = J £,{r)e lk - r d 3 r (70) 

and with inverse transform: £(r) = (2ir)~ 3 J P(k)e~' lk ' r d 3 k. At the origin r = 
we obtain 

&(Q) = (&(*)) = J P(k) W 2 (kR)d 3 k (71) 

where W is the window function function that reflects the filtering of the field. 
The power spectrum is the contribution of modes of wavenumber k to the total 
variance, per unit volume of wavenumber space. If the fluctuation field is a 
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Gaussian Random Field, then the power-spectrum contains all the statistical 
information of the fluctuations. A similar formulation is applicable also in the 
two dimensional case (where the density field is on the surface of a sphere - the 
sky) . Only that instead of the Fourier transform we use the Spherical Harmonic 
transform. 

Although I will not enter in the details of how to estimate the power spec- 
trum of some distribution of extragalactic objects, I will only note that a good 
estimation of the window function (containing the survey boundaries, obscura- 
tion, radial selection function and instrumental biases for example) is necessary 
in order to get a reliable power-spectrum determination. 

A further ingredient that complicates considerably matters is the so called 
biasing of the galaxies, or in general of any mass-tracer population, with respect 



to mass 1 79 . The usual relation, assumed between the mass-tracer fluctuations 
{Sti) and the underline mass fluctuation field, is encapsulated in the bias factor 
b: 

Str = b 5 m (72) 
and therefore we have that the galaxy (tracer) power spectrum is 

P tr (fc) = b 2 P m (k) 

It has been shown that the linear biasing model (|T2|) is a good approximation, 



at least on scales were non- linear gravitational effects are weak (see [106| and 
references therein). 

Within the inflationary paradigm the initial fluctuations, in the early uni- 
verse, that gave rise to the observed large-scale structure today, are adiabatic 
and Gaussian and therefore one can characterize these fluctuations completely 
using the above tools. The power spectrum of such initial fluctuations is: 

P in (k) = Ak n 

usually with n = 1 (Harris on- Zeldovich spectrum), and A its amplitude. The 
different fluctuation damping mechanisms, operating during the radiation dom- 
inated area, modify P; n (fc). These effects can be encapsulated in the transfer 
function, T(k) and today's linear fluctuation spectrum has the form ||: 

P(k)=T 2 (k)P in (k) 

In the linear regime (while fluctuations 1) the power spectrum shape is pre- 
served, because each Fourier mode evolves independently. 

I will now concentrate on a few methods, based on the clustering of galaxies 
and QSO's, and I will present only very recent results. The subject is extremely 
rich, many have laboriously worked towards attaining the goal of pining down 
the different cosmological parameters and I hope that they will forgive me for 
not being able to mention the vast literature on the subject. 
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Fig. 10. Pie-diagram of the latest 2dF release, containing 160000 galaxies (with per- 
mission from the 2dfgrs team). 



The Shape of P(fc): The popular Cold Dark Matter (CDM) model has a T(k) 
parametrised by the so-called shape-parameter, r, which characterizes the shape 
of the P(k) and has the form [168]: 



r = Q m h exp -J? B (1 + V2h/fin 



(73) 



We see that measuring the the power-spectrum of extragalactic populations and 
estimating r we can put constraints on the combination of the cosmological 
parameters: fi m and h (Qb affects weakly T). 

The recent 2dF survey |J7| has measured already measured more than 160000 
galaxy and 10000 QSO redshifts, which constitute it the largest spectroscopic 
catalogue of extragalactic objects. A number of recent papers have estimated 
the galaxy and QSO power-spectra providing important constraints on the r 
parameter. From the P(k) of galaxies [ 1 1 7 1 it was found that there is a degeneracy 
between the f2 m h and 12b //2 m , which if broken (using the CMB or BBN results) 
provide the following constraints: 

Q m h ~ 0.2(±0.03) f2 B ^ 0.15(±0.07)J? ro 

which for h = 0.72 means that fl m = 0.28. The analysis of the QSO P(k) 
[ [70| showed a somewhat smaller value but still consistent within the statistical 
uncertainties: 

r~o.i(±o.i) 
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Similar results come from the SSDS photometric survey |165] which contains 
1.5xl0 6 galaxies with redshifts up to ~ 0.4. The angular power-spectrum anal- 



ysis, after inverting to 3D using Limber's integral equation, give |47j] [17C] 



r~o.u±° ;ll 95% ex. 

Obviously, all these results support a low J7 m universe. 

Redshift Space Distortions: As we have already discussed in section 2, the 
measured expansion velocity of an extragalactic object contains also the con- 
tribution of the local gravitational field. The peculiar velocities will distort the 
apparent 3D distribution of the extragalactic objects, a fact which will mani- 
fest itself in the 2-point correlation function, when plotted as a function of the 
transverse and radial pair separation. The redshift space correlation function is 
related to the real space, under a few assumptions; see J78|, according to: 

&M = 6M + + (74) 

where (3 = fi^f /b. By estimating the angular power-spectrum, Ci, or its Fourier 
transform w(6), and then inverting it to 3D via Limber's equation, one has an 
unaffected, by redshift space distortions, measure of these parameters. Then 
using (|74|) and the measured £ s (r) one can place constraints on j3. Such an 
analysis of the 2dF galaxy survey gave |113| : 

Q m ~ 0.25(±0.06) 6- p 166 

A subsequent analysis, using a different method, provided very similar results 
III: 

Q m ~ 0.3(±0.15) 6-p 1 ' 66 



The corresponding QSO survey [ 108 1 did not provide very stringent constraints 
(ic., the EdS model was rejected only at a 1.4ct level), however their best fit 
gives: 



1.66 



n m ~ 0.21(±0.15) b Q ^ 

Joint Likelihoods: Alot of recent interest was generated by the understanding 
that joining the analyses of different data sets, one may break the degeneracies 



between the different cosmological parameters (cf. [182] [[33| |29J |181|). Espe- 
cially joining the CMB, SNIa and large-scale clustering results may lead to strong 
constraints on more than 8 cosmological parameters (see however [^7| for many 
subtleties involved). 

The joint analysis [[34] of the 2dF galaxy P(k) and the CMB data have pro- 
vided another strong indication for a positive cosmological constant, independent 
of the SNIa results, with a 2a range: 

0.65 < i? A < 0.85 

Furthermore, some of the other constraints are: 0.17 < J? m < 0.31 for h = 0.72, 
0.6 < h< 0.86. 
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3.4 M/L observations 

Each astronomical object is characterized by a quantity called Mass-to-Light 
ratio, M/L. A convenient scaling of M/L is done by using the value of the 
solar neighbourhood, Mq/Lq. Then, observing M/L > Mq/L q would imply 
the existence of DM of unknown composition and origin. Note that in most 
cases the evidence for M/L > Mq/Lq comes from the gravitational effects of 
the DM. Different classes of extragalactic objects (galaxies, clusters, etc) are 
characterized by different M/L, indicating that possibly a different fraction of 
the total mass of each type of object is Dark. 

Estimating the universal luminosity density (see further below) and using the 
derived M/L values of each class of extragalactic object we can estimate its con- 
tribution to the total i? m . Furthermore, if the estimated M/L is representative 
of the global universal value, then we can derive the overall value of J? m . 

I will present below the basic ideas behind the determination of M/L for the 
different extragalactic populations. 

Spiral Galaxies: The Rotation curves of spirals (see Fig. 10) are obtained by 
measuring Doppler-shifts of emission lines in HII regions, at radio wavelengths 
using the 21-cm emission line of neutral Hydrogen or using the CO-line and in 
the latter cases the rotation curve is mea sured at several times the optical radius 
of a galaxy (for a recent review see fl62|| ). 

One would expect the rotation curve to fall roughly as the surface brightness 
and beyond a few length-scales to fall as v mt oc r -1 / 2 (because most of the mass 
is rather centrally located) which is not observed. The rotation curves are found 
to be flat as far as they can be observed. From simple Newtonian Physics we 
have that: 



and since v mt oc constant, we have that M(r) oc r, ie., mass increases linearly 
with distance beyond its optical radius, an indication for the presence of dark 
matter. The average value of M/L found for spiral galaxies out to ~ 20 kpc, is: 



Note that the M/L value is an increasing function of outer radius, implying the 
existence of an extended dark matter halo (cf. ||). 

Elliptical Galaxies: In principle one can invoke a similar method with that 
of the spirals, if instead of the rotational velocity, the stellar velocity dispersion 
is used. The amplitude of this velocity measure is dictated by the gravitational 
potential of the elliptical and thus the virial theorem can be used to determine 
its total gravitating mass. However, stellar orbits are highly uncertain and the 
velocity dispersion ellipsoid may not be isotropic (see discussion in |jl6[ ). 

Luckily, elliptical galaxies contain very hot gas (T~5x 10 6 K) which emits 
X-rays. This gas constitutes about 10% of the observed mass in stars. For a 




M/L ~ 10 ± 2M /L© 



(75) 
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Fig. 11. R otation curves of 4 spiral galaxies, one of which is the Galaxy (data taken 



from [161 



spherically symmetric galaxy in hydrostatic equilibrium (v = 0) we have from 
Euler's equation that: 

+ • V) v = — Vp - V<5 =$> 
at p 



dp 
dr 



GM(< r)p 



Using the ideal gas law (j> = pk^T/m) with m the molecular mass and fcs the 
Boltzmann constant, we obtain that: 



M(< r) = 



k^Tr f dlnp dlnT 



pm p G 



d In r d In r 



(76) 



where T is the gas temperature, m p is the proton mass and p is the mean 
molecular weight. Therefore if we measure the temperature and density profiles 
we can find the total mass distribution M(< r). Finally the average mass-to-light 
ratio from this class of objects, out to ~ 20 kpc, is: 



M/L ~ 25 ± 5M Q /L Q 



(77) 



Groups of Galaxies: Groups of galaxies containing a few galaxies (3 - 10) are 
usually considered bound (due to the high density of galaxies) . Then, according 
to the Virial theorem (2T + U = 0), a group with TV galaxies having velocity 
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and position vectors (relative to the centre of mass of the group) v\ and r\ 
respectively, has: 

N Gmm- 

i=l i=l j<i 1 1 Jl 

However what we measure is not Vi but only the line-of-sight component of the 
the velocity u\. Assuming isotropic orbits then Assuming that the 

mass-to-light ratio of each member galaxy is the same, mj/ Li = M / L, we obtain: 

M = Sir (m 
L 2GEE j <i^j| J Ri-«jh 1 1 1 

where \R[ — Rj \ is the projection of \ri — rj| on the plane of the sky. 

However this measure should be used only if N is large, or otherwise sta- 
tistically, ie, averaged over a large number of groups to find the (M/L) of the 
ensemble. Note that this estimator is very sensitive to the group-finding algo- 
rithm, selection procedure, selection function corrections (ie., corrections to take 
into account the fact that the galaxy density artificially decreases with depth in 
magnitude/flux limited samples). Furthermore, projection effects (ie; inclusion 
in the group of a foreground or background galaxy) is also a source of error. 
For example the inclusion of a non-member would result in an artificially higher 
M/L while conversely the exclusion of a member (due to a high iti) would un- 
derestimate the M/L. 

The mass-to-light ratio for groups of galaxies (cf. |l32fl ) is: 

M/L ~ 180 ± 50M Q /L Q . (79) 



Clusters of Galaxies: The most common approach to measure the cluster 
M/L is based on the assumption that clusters are in virial equilibrium, for which 
we do have strong indications. The first step in determining M/L is to measure 
the cluster velocity dispersion accurately. Using the recent compilation of 395 



R > ABELl/aco cluster velocity dispersions [167] I find 



cr = ((« - vf) 1/2 S3 750 ± 322 km/sec 

However, this value is based even on clusters for which the velocity dispersion 
was determined from a very small number of galaxy redshifts, a fact which 
increases the possibility of assigning an erroneous cr- value. This can be seen in 
the left panel of Fig. 11 were we plot the normalized a frequency distribution for 
clusters with N z < 20 and N z > 20. The skewed distribution to large cr-values 
for N z < 20 is evident, while the mean cr- value is smaller than for the N z > 20 
case. As a compromise between having enough redshifts per cluster, in order to 
get a reliable cr-value, and a large cluster sample, we choose those with N z > 20. 
For these 195 clusters I find that: 



a w 822 ± 294 km/sec 
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Fig. 12. Left panel: Frequency distribution of ABELl/aco cluster velocity dispersions 
based on N z < 20 (empty histogram) and on N z > 20 (shaded histogram) . Right panel: 
The corresponding M/L distribution. 



A consistent value of a ~ 940 ± 208 km/sec has been found from the 16 high 



redshift clusters of the CNOC project |32| and a variety of methods used to 
define cluster membership and account for interlopers p3| . 
Now from the Virial theorem, M = 2a 2 r a /G, we have: 

M c (< r a ) « 7 x 10 14 (^) 2 h- 1 M Q . (80) 

where r a = 1.5 hr 1 Mpc. Using a number-galaxy weighted luminosity estimation, 
and assuming an average cluster value in the optical of L ~ 10 12 hr 2 L Q (cf. 
[|J), we obtain the distribution of M/L values (see right panel of Fig. 11). Since 
the distribution is non-Gaussian we quote below the median and 68% confidence 
levels: 

(M/L) c w 3201^ h Mq/Lq (81) 

A more accurate method that can be used to estimate cluster mass is based on 
the measurements of the X-ray emission from the ICM gas, similar to the method 
used for individual elliptical galaxies ((76]). Although not strictly correct, usually 
an isothermal cluster profile is used (dlnT/dlnr = 0), which greatly simplifies 
calculations. However, recent experiments have shown that indeed the cluster 
Temperature does not vary significantly with radius (dlnT/dlnr ~ —0.7 ~ 
—0.8). Estimates of M/L derived with this method are in general agreement 
with 
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Table 2. Mass-to-light ratios and contribution to J? m for the different scales. 



Object 


scale (h' 1 Mpc) 


(M/L) hr L 




Spirals 


0.02 


10 ±2 


0.0071 ±0.0015 


Ellipticals 


0.02 


25 ±5 


0.018±0.004 


Galaxy pairs 


O.f 


80 ±20 


0.057 ±0.012 


Groups 


0.8 


180 ± 60 


0.13 ±0.09 


Clusters 


1.5 


3201J™ 


23+ 12 



Global luminosity density and f2 m : We estimate the value of (L) using the 
galaxy luminosity function which is defined such that <P(L)dL is the number 
density of galaxies having total luminosity in the interval (L, L + dL). A good fit 
to the observed luminosity func tion of the field population of galaxies is provided 
by the Schechter function |l53| : 



W = ^(£) Q exp[-L/L.] (82) 

where a = —1.29 ± 0.11, (f>* = 1.3 ± 0.3 x lO^h 3 Mpc" 3 and L* = 1.1 x 
10 10 h~ 2 L Q (see also (52)). Then the mean luminosity density, corresponding to 
13), is: 

(L) = J L$(L)dL = ^L„r(a ± 2) w 2 x l0 s hL Q Mpc" 3 (83) 



1400 n m h Mq/L g . (84) 



Using the value of p cr (|15|) we have: 

M p ftmPc 



L (L) (L) 



In Table 2 we summarize the (M/L) values found at different scales and the 
corresponding contribution to i7 m . 

Since galaxy clusters are the deepest potential wells in the Universe and they 
accumulate baryonic and DM from large volumes, it is expected that their M/L 
ratio could represent the Universal value. This view is supported by the fact 
that the increasing trend of M/L with scale (seen in Table 2) reaches a plateau 
at the corresponding value of the clusters (cf. ||, Q). Therefore the universal 
value, as given by the clusters is: 

(7 ~ oq+ 012 

m — u ° — 0.06 

The analysis of the CNOC sample of 16 distant clusters (0.17 < z < 0.55), 
provides a consistent value but with significantly smaller uncertainty; f2 m ~ 



0.19 ±0.06 (see [32 and references therein). 



3.5 Cluster Baryon Fraction: 

Assuming that on the scales from which clusters accrete matter during their 
formation, there is no segregation of baryons and DM, then the ratio of baryonic 
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to total matter in clusters is representative of the universal value ((2b /^m), 
something supported also by hydro-dynamical numerical simulations. Galaxies 
add up to only ~ 5% of the total cluster mass, while the hot gas, which fills the 
space between galaxies, accounts for ~ 20%. 

So adding up the galactic and gas contribution to the mass of the cluster 
we obtain a measure of the total baryonic cluster mass, Mb and then if we can 
measure the total mass of the cluster, M tot , we can estimate i? m . 

Assuming hydrostatic equilibrium, the gas tra ces the clu ster total mass, and 



using ([76D we can obtain, M to . Most studies (cf. |86| [103| |5J]) find: 

Me 



0.10-0.13 for h = 0.72 (85) 



From the primordial nucleosynthesis constraints (see p3| ) and from the recent 
Bommerang results we have that H B h 2 ~ 0.020 ± 0.004), which for h = 0.72 
gives J?b — 0.04 and therefore: 

Q m ~ 0.35 (±0.05) 



Note however that in a recent study [147|, the application of various corrections 



to account for the clumping of gas and the gas fraction gradients, within the 
virial radius of a cluster, resulted in significantly lower values of the baryon 
fraction and thus higher values of J? m . 



3.6 Large-Scale Velocity Field: 

As already discussed in section 2, the local gravitational field produces peculiar 



velocities superimposed on the general expansion (35). Measurements of peculiar 
velocities can provide direct information on the mass content of the Universe, 
since they can be related to the density fluctuation field, which itself can be 
observed directly. The basic idea is that a body of mass M will produce a different 
gravitational field if it is embedded in a low or high density Universe. In a high 
density Universe it will correspond to a lower density fluctuation than in a low 
density Universe and thus it will produce weaker/stronger gravitational effects, 
respectively. 

If S(x) — (p(x) — p)/p is the mass fluctuation at x then using linear pertur- 
bation theory, continuity, Euler and Poisson equations, we obtain the that the 
growing mode of the evolution of fluctuations is: 



5 (x D{t) 



5 _D 
5~D 



Note that in the EdS universe we have D(t) — 1 3 . Linearizing the mass continuity 
equation, we obtain: 

Vv = -R5 = R8— (86) 
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which has solution: 

where the the dipole vector D(x) is related to the gravitational acceleration 
vector and b is the bias factor relating the mass-tracer fluctuations (<5tr) with 
the underline mass fluctuation field (|72|). Note that in most cases it is expected 
that b > 1, ie., the tracer distribution is more clustered than the matter. This has 
been born out from the study of random Gaussian fields in which the higher the 
density fluctuations the more clustered they are (cf. § ) and from the correlation 
function analysis of extragalactic objects by which it was found that the relative 
bias factor of ABE LL clusters, optical and IRAS galaxies is [b c \ : 6 op : &ir = 4.5 : 



1.3 : 1] (sec [ 111 ). ft is therefore expected that this hierarchy of decreasing 



correlations shoul d con tinue to the underlying mass distribution. Note that (| 



can be written as 114 



" = 5W < 88 > 

One then needs good estimates of peculiar velocities, knowledge of the matter 
tracer fluctuation field to estimate D(x), and an understanding of the biasing 
between matter and light in order to put constraints on J? m . 



Local Group Dipole: The above test was first applied to the Local Group, 
since from the CMB dipole we have an excellent measurement of its peculiar 
velocity. The dipole moment, D(x), is measured using different populations of 
extragalactic objects (iRAS, optical galaxies, AGN's, ABELl/aco optical or X- 
ray clusters), weighting each object by r~ 2 or by a weight proportional to r~ 2 
(like flux, diameter 2 ). Under the assumption that light traces the mass, then the 
dipole moment is a measure of the peculiar force acting on the LG. 

In linear theory the peculiar velocity is parallel to the acceleration and there- 
fore finding an approximate alignment of the two vectors tells us that the fluctu- 
ations causing this motion are present within the depths of the sample, provided 
that the dipole moment converges to its final value before the characteristic depth 
of the sample. Furthermore, it tells us that possible local non-linear effects do 
not strongly affect the D(x) determination. 

The multipolc components of some mass-tracer distribution are calculated 
by summing moments. For example the monopolc and dipole terms are: M — 
l/47r^Wi and D = 3/47r^Wiri where Wi cx r~ 2 . The dipole vector, D, is cal- 
culated by weighting the unit directional vector pointing to the position of each 
galaxy, with the weight Wi, of that galaxy and summing over all N available 
galaxies in the survey. Note however, that even a uniform distribution would 
produce a dipole, the so-called shot-noise dipole, if it is sparsely sampled. The 
shot noise error is (D-D) 5 = 3/4tt N2(w 2 )2, and even in a clustered dis- 
tribution, the estimated dipole will always have such a shot-noise contribution 
which should be taken into account. A further difficulty is that whole-sky dis- 
tributions of extragalactic mass-tracers are unavailable, either due to survey 
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Fig. 13. Dipole amplitude build up as a function of distance of the IRAS PSCZ galaxy 
and the ABELl/aco cluster samples. The ABELl/aco cluster dipole has b een scaled 
down by a factor of ~ 4, to take into account the relative bias factor (from [124]). 



limitations, extinction near the Galactic plane, cirrus emission, magnitude, flux 
or diameter limits (which will cause a different population mix to be sampled 
at different depths, and the contributions of the 'faint' or 'small' objects to 
be missed). Therefore, complicated corrections should be applied, which usually 
have well understood properties and in any case are always tested with numerical 
or Monte-Carlo simulations. 

Results from many different galaxy catalogues (optical or IR) tracing depths 
up to ~ 100 — 120 h~ 1 Mpc, show dipoles that are well aligned with the CMB 
dipole (mass dipole) which indicates that galaxies trace the mass distribution. 
Some recent analysis of IRAS galaxies provide values: 

Q m « 0.62(±0.13) b\i G 



(cf. [166], [145 and references therein) and 



n m « 0.55(±0.2) blf 

(cf. [Q). The difference between the above values is most probably due to the 
different biasing factors ie., bm < 6 op , which is natural since optical galaxies 
trace more accurately the deep potential wells (clusters) while the IRAS galaxies 



trace better the field, fn fact it has been found that b op /bm ~ 1.2 — 1.4 (cf. [ 1 1 1 1 
@ jl|§ and refer ences therein). Furthermore, if IR galaxies are biased with 



respect to mass (ie., if &ir > 1.2) then the above results could imply: J? m ~ 1. 
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Now, galaxy clusters, being the largest gravitationally-collapsed structures in 
the universe and luminous enough to be detected to very large distances, have 
also been used to probe the local acceleration field. Most studies are based on the 
optically selected ABELL / ACO clusters and they provide strong evidence that 
the LG dipole has significant contributions from depths up to ~ 160/i _1 Mpc 



(see Fig. 12 and [152], 122], pa]). However, due to the the volume incomplete- 



ness of richness class R=0 clusters (cf. [ 1 1C ] ) and to optical projection effects 
(enhancement of galaxy density along the direction of foreground rich clusters 
which cause inherently poor background clusters or groups to appear rich enough 
to be included in the sample), these results had to be checked. Using well de- 
fined X- ray c luster samples free of the above effects [[[l| indeed the results were 



verified 123 . The results based on optical or X-ray clusters, imply: 



il m « 0.07 - 0.09 bif 6 



Note that from the correlation function analysis of ABELl/aco clusters [111] we 
have that the relative bias factor between clusters and IRAS galaxies is around 
6 c ijr, ~ 4.5 and if the bias factor of IRAS galaxies with respect to the mass is 
b g > 1, then the above results are consistent with Q m ~ 1. 

It is extremely interesting that the galaxy distribution in the local universe 
(within ~ 100 hr 1 Mpc) produces a gravitational acceleration that is aligned 
with that produced by the matter distribution (as determined by the CMB 
dipole) and in the same time the distribution of clusters of galaxies on larger 
scales R < 250 h~ x Mpc produce also a gravitational acceleration aligned with 
that of the matter distribution. Furthermore the galaxy and cluster distributions 
produce Local Group acceleration profiles that are directly proportional to each 
other (see Fig. 13), while the dipole in equal volume shells, seems to be roughly 
aligned with the CMB dipole direction out to very large depths (cf. Jl22| , jll| ) . 
This implies that there is a coherent anisotropy in the mass distribution over a 
region with a diameter of ~ 300 Mpc, which then sets a lower limit to the 
scales over which the Cosmological Principle applies. 

However, there is a dichotomy among different studies trying to identify the 
convergence scale of the dipole, or equivalently what is the largest scale over 
which we observe bulk motions. The previously discussed dipole studies as well 



as some peculiar velocity studies (cf . [|39j , JT5J , ||188| ) , support the view of a large 
convergence depth (radius of ~ 150 h^ 1 Mpc). Other peculiar velocity studies 
(cf. |3q] , |Il| |, ]40|] , |3t| , ]39[| ) support a significantly smaller convergence depth 
- 60 h- 1 Mpc. 



POTENT - from Radial Velocities to Density Field: In the previous 
analysis only one velocity was used, that of the Local Group, and although it 
is very well measured we still have the problem of cosmic variance. Therefore, 
ideally the velocity-acceleration comparison should be performed for a number 
of "observers". This has been possible due to the POTENT algorithm proposed 
in |l4| and developed extensively by Dekel and his collaborators (cf. jl6) and 
references therein). 
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Fig. 14. (Included separately as a JPG image) Comparison of predicted density and 
velocity fields: Left panel shows the observed ABELl/aco density field and the corre- 
sponding predicted velocity field. Right panel shows the matter density field predicted 
by POTENT and the peculiar velocity field of Mark III galaxies (from Q). 



The basic idea follows. The large-scale velocity field, evolving via gravita- 
tional instability, is expected to be irrotational Vxu = 0. This remains a good 
approximation in the mildly-nonlinear regime as long as the field is properly 
smoothed. This implies that the velocity field can be derived from a scalar po- 
tential, 

v(x) = -V$(x) , 

and thus the potential can in principle be computed by integration along the 
lines of sight, 

&(x) = - I u(r',9,(j))dr' . (89) 
Jo 

The two missing transverse velocity components are then recovered by differen- 
tiation. Then from ( |S6| ) we recover the density fluctuation field, which can then 
be compared to the observed density field, determined from large whole-sky sur- 
veys. The current sampling of galaxies enables reliable dynamical analysis, with 
a smoothing radius as small as ^10/i _1 Mpc, where |V-t>| obtains values larger 
than unity and therefore mildly non-linear effects play some role. 

The most reliable density-density analysis, incorporating certain mildly non- 
linear corrections, is the comparison of the IRAS 1.2 Jy redshift survey and the 
Mark III catalogue of peculiar velocities yielding, at Gaussian smoothing of 



12/i- x Mpc 1 15? 



1.66 



Q m ~ 0.82(±0.16) b\£ 

A similar analysis, using optical galaxies J74| has provided a somewhat lower 
value: 

n m ~ 0.6 (±0.15) bl f 

in accordance with the expected higher biasing parameter of optical galaxies 
with respect to IRAS ones. These results are consistent with the dipole analyses 
and with J? m ~ 1 for bm > 1.2. 



However, a variety of methods using v-v comparisons (eg. VELMOD - [187]) 



developed to compare observed and derived velocities (using either the IRAS or 



ORS gravity fields), typically yield values of: (cf. [189] and references therein, 



n m ~ 0.3 (±o.i) &! 1 - 66 



Q m ~ 0.14 (±0.05) bl f 

which are consistent with J? m < 1 for any reasonable value of 6 g . Therefore, there 
seems to be a discrepancy between different analyses, even if in some cases, they 
use the same data, a fact that needs further study and tests of the reliability of 
each method. 
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A study [^7j using the ABELl/aco clusters to trace the density field and 
comparing it with the potent reconstructed field from the Mark III catalogue 
of peculiar velocities (see Fig. 14) found: 

Q m « 0.07 - 0.09 blf 6 

in good agreement with the dipole analysis of ABELL / ACO clusters and consistent 
with f2 m ~ 1 for the estimated value of b c \ [111]. 



Local Group infall to Virgo: This is an interesting method to calculate i? m 
on scales of ~ 10 h^ 1 Mpc. One relates the Local Group infall (towards the 
centre of the Local supercluster) velocity with the acceleration induced to the 
LG by the mass overdensity in the Local Supercluster, assuming a point-mass 
approximation. We have from (87): 



2 GSM 



3 ffoffl 4 



(90) 



From (|72j) we have that the galaxy fluctuations is a biased tracer of the underline 
mass fluctuation field: 



SM 



5M 



M5„ 



(91) 



Thus from (90) we have: 



i a - 1 



3 b 



cz S g (l + 6 g r^ 



where the mildly non-linear correction on the right-hand side is according to 
JL Q 1| . A recent study [174| using the SBF method to determine the local velocity 
field within cz < 3000 km/sec find a Virgo-centric infall of v- m ~ 140 km/sec in 
agreement with ~ 160 km/sec, implied from the Virgo contribution to the X-ray 
cluster dipole [ 123 1 . Furthermore we have that cz ~ 1005 km/sec |104[ Jl49| and 
<L ~ 2.8 ± 0.5. Therefore we obtain: 



n m ~ o.i(±o.o2) 



150 km/sec 



1.66 



a.66 

op 



Although this method of determining f? m is 'clean', the fact that the local pe- 
culiar velocity field is affected by mass concentrations well beyond the Local 
supercluster introduces a further uncertainty in the determination of f2 nl . 



Velocity-Field results: The outcome of the different large-scale dynamical 
studies do not converge to a unique value of the mass density parameter. There 
is need to check the methods and understand the source of this discrepancy. In 
Table 3 I sum up the different results from the different velocity-field analyses. 
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Table 3. Results from some recent velocity field studies 

Type of study fl m y M fl m ^ p LM ^ m b^ 1 

Dipole 0.62 0.55 0.08 

POTENT/Mark III 0.82 0.6 0.08 

VELMOD 0.3 0.14 
LG-infall 0.1 



3.7 Rate of Cluster Formation Evolution: 

The rate of growth of perturbations is different in universes with different matter 
content. For example, the perturbation growth in a 4? m = 1 universe is propor- 
tional to the scale factor, ie., 6 oc (l + z) -1 , while in the extreme case of an empty 
universe (Q m = 0); 8 = constant. From ( |32| ) we see that J? m < 1 universes will 
behave dynamically as an fi = 1 universe at large enough redshift, and at some 
redshift z ~ 1 curvature dominates and perturbations stop evolving and freeze, 
allowing clusters to relax up to the present epoch much more than in an J2 m = 1 
model, in which clusters are still forming. This can be seen clearly in Fig. 15 
where we plot the evolution of the perturbation growth factor, defined as: 

din 5 
1 dlni? ■ 

For a n A = universe, / ~ (cf. jUj ) and thus for the EdS model / = 1. For 
J7a > models there is a redshift dependence of /, but in the present epoch it is 
indistinguishable from the corresponding value of the open (f2 m = 1 — J?a) model 
[ |S6| . It is evident that an < J7a < 1 universe behaves as an J? m = 1 model up 
to a lower redshift than the corresponding open model, while at redshifts z < 1 it 
behaves like an open model, which implies that clusters should be dynamically 
older in such a model than in the EdS. Therefore one should be able to put 
constraints on fl from the evolution of various indicators of cluster formation, 
especially in the range where the dynamical evolution between the models differs 
maximally (vertical dashed lines in Fig. 15). 

Ideally, one would like to study the evolution of the cluster mass function 
but since light is what we observe (temperature as well - due to the hot ICM 
X-ray emission) , various related indicators are usually studied (Luminosity func- 
tion, temperature function, morphology etc), but then one has to pass through 



the machinery provided by the Press-Schechter formalism [128 , which gives the 
mass function of collapsed halos at any epoch as a function of the cosmological 
parameters that enter through the assumed power spectrum of perturbations. 

Luminosity function: Based mostly on Einstein and Rosat surveys, many 
studies have found an evolving X-ray luminosity function, ie., less z > 0.3 clusters 
than expected for a no-evolving luminosity function, ie., a negative evolution (cf. 
[ |S3"| and references therein). Such a behaviour is expected in models with 

n m ~ i 

However, see [|9| for a different view. 
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Fig. 15. The evolution of the perturbation growth factor, /, in 3 models (EdS, open 
i? m = 0.3 and the currently popular flat Ha = 0.7). The vertical dashed lines indicate 
the redshift range 0.3 < z < f where the 3 models, jointly, differ maximally. 



Temperature function: Estimates of the temperature of the X-ray emitting 
ICM gas can be reliably estimated from the iron line-emission. Then the cluster 
temperature can be either transformed to a mass (assuming hydrostatic equi- 
librium and isothermality) and thus deri ve a mass function to compare with 
the Prcss-Schechter predictions (cf. 133 and references therein) pointing to 
/2 m < 0.3 , or use the evolution of the temperature distribution function. Again 
different studies find either no evolution (cf. 0]) pointing to 



On 



0.3-0.5 



or evidence for evolution 18C] |f9(] pointing to 

i? m -0.7-1 



Evolution of L — T relation: Under the assumption of hydrostatic equilib- 
rium and isothermality one can easily show, from Euler's equation, that the 
bremsstrahlung radiation temperature is T oc M v / R v (where M v and R v are 
the c luster virial mass and radius). Using the spherical collapse top- hat model 
[ pT4| one obtains R v oc T 1 / 2 A(z)- 1 ? 2 E(z)- X l 2 , and then by using @: 

L x ex / g 2 Af v 2 T 1 /2 i? -3 ( 92 ) 
where f s is the gas mass fraction. Then one finds (cf. pf[|): 



L x oc f z T 2 A{z)V* E{zfl 2 
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where E{z) is given by (|l^) and A(z) is the ratio of the average density within the 
virialized cluster (< R v ) and the critical density at redshift z, which also depends 
on the cosmological model. However, this model fails to account for observations 
which show a steeper T-dependence, L x oc T 3 " 3 3 (such a dependence can be 
recovered from (|9^) if f g oc T 1 / 2 ). In any case, the L — T relation is expected 
to evolve with time in a model dependent way. Most studies (see references in 
[ p.55| ) have found no evolution of the relation while a recent study of a deep 
(z ~ 0.85) Rosat cluster survey |2fJ found: 

Q m ~ 0.35 

However, there are many physical mechanisms that affect this relation (eg. gas 
cooling, supernova feedback etc) and in ways which are not fully understood (cf. 

H). 



Evolution of Cluster Morphology: As we have already discussed, in an open 
or a flat with vacuum-energy contribution universe it is expected that clusters 
should appear more relaxed with weak or no indications of substructure. Instead, 
in a critical density model, such systems continue to form even today and should 
appear to be dynamically active (cf. 134], |59| , f85| ). Using the above theoretical 
expectations as a cosmological tool is hampered by two facts (a) Ambiguity in 
identifying cluster substructure (due to projection effects) and (b) Post-merging 
relaxation time uncertainty (cf. [ 15C ] ) . However, criteria of recent merging could 
be used to identify the rate of cluster morphology evolution and thus put con- 
straints on O m . Such criteria have been born out of numerical simulations (cf. 
[ fL39| , [140 1) and are based on the use of multiwavclcngth data, especially optical 
and X-ray data but radio as well (cf. [192], [ 1 54 1 ) . The criteria are based on 
the fact that gas is collisional while galaxies are not and therefore during the 
merger of two clumps, containing galaxies and gas, we expect: (1) a difference 
in the spatial positions of the highest peak in the galaxy and gas distribution, 
(2) due to compression, the X-ray emitting gas to be elongated perpendicularly 
to the merging direction, and (5) temperature gradients to develop due to the 
compression and subsequent shock heating of the gas. The first two indicators 
are expected to decay within ~ 1 Gyr after the merger, while the last may sur- 
vive for a considerably longer period (see for example Fig. 16) . For such a study 
to be fruitful, a large number of clusters, ideally covering the redshift range 
0-3 J$ x < 1, must be imaged in both the optical and X-ray band. 

However, a rather cruder but still useful test of cluster morphological evo- 
lution could be used. For example, cluster ellipticity is a relatively well defined 
quantity; although systematic effects due to projections in the optical or the 
strong central concentration of the X-ray emitting gas (since L x tx n 2 ), s hould 
be taken into account (cf. [pTJ ) . An early study, using the Lick map 121 , had 



found that cluster ellipticity decreases with redshift, however due to possible sys- 
tematic effects involved in the construction of the data, they did not attach any 
weight to this discovery. Recently, two studies using optical and/or X-ray data 
[ |100| [127 1 (see also [ 120 1 ) found that indeed the cluster ellipticity decreases with 
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Fig. 16. Optical apm (colour) and Rosat X-ray (contour) images of 2 abell clusters. 
Peaks of the APM galaxy distribution is shown in blue. A3128 has the signature of a 
recent merger: the peaks in the distribution of galaxies and in X-ray emitting gas are 
orthogonal to each other. A2580 on the other side seems a sm ooth relaxed cluster with 
the gas and galaxies tracing the cluster potential (from [ 12C ]). 



redshift in the recent past, z < 0.15 (see Fig. 17) This was interpreted by [10C] as 
an indication of a low-/2 m universe because in such a universe one expects that 
merging and anisotropic accretion of matter along filaments will have stopped 
long ago. Thus the clusters should be relatively isolated and gravitational relax- 
ation will tend to isotropize the clusters reducing their ellipticity, more so in the 
recent times. 

If this is the case then one should expect an evolution of the temperature of 
the X-ray emitting gas as well as the X-ray cluster luminosity which should fol- 
low the same trend as the cluster ellipticity, decreasing at recent times, since the 
violent merging events, at relatively higher redshifts, will compress and shock 
heat the diffuse ICM gas [ 138 1 . Such evidence was presented in 127 using a com- 
pilation of measured ICM temperatures and luminosities in two volume limited 
X-ray cluster samples (based on the XBAC and BCS samples). Also, one could 
naively expect an evolution of the cluster velocity dispersion, increasing at lower 
redshifts, since virialization will tend to increase the cluster 'thermal' velocity 
dispersion. In |55) no evolution was found between a local sample (z < 0.15) and 
a distant one 0.15 < z < 0.9. However, unrelaxed clusters can also show up as 
having a high velocity dispersion due to either possible large peculiar velocities of 
the different sub-clumps [142] or due to the possible sub-clump virialized nature. 



Therefore, a better physical understanding of the merging history of clusters is 
necessary in order to be able to utilize the velocity dispersion measure as an 
evolution criterion. 

Other related studies, using the morphological characteristics of the large- 
scale structures, have been used to place cosmological constraints. For example, 
the shapes of superclusters and voids, using the IRAS-PSCZ redshift survey and 
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Fig. 17. The evolution of ellipticity in APM clusters with significant substructure [120 1 



the ABELl/aco cluster distribution show a clear preference for a A-CDM model 
over a Q m = 1 model @ |2§ f| (l02) but see @. 
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4 Summary 

I have attempted to present the basic ideas behind many (but not all) of the 
current methods used to estimate various Cosmological parameters, especially 
the Hubble constant (H ), the curvature of the Universe {Q) and its matter 
content (i? m ). To do so I reviewed some basics of the standard, Robertson- 
Walker, Cosmology in order to highlight the interrelations of the Cosmological 
parameters and the way they affect the global dynamics of the Universe. 

The results of a variety of different analyses, based on a multitude of data, 
point towards a concordance model, which is a flat, £2\ ~ 0.7, with an inflationary 
spectral index n — 1, h ~ 0.72, J?b — 0.04, t ~ 43. f Gyr's. Nevertheless there 
are conceptual problems and open issues that may or may not prove to be 
daunting. 

For example, the J? a = 0.7, h — 0.7 model seems to have problems generating 
the correct power on galactic scales. Detailed numerical simulations of [|o) show 
that this model has much more power on small scales to be reconciled with obser- 



vations (see discussion in [12E]). Furthermore, there is the fine tuning problem. 
Why does the energy density of the vacuum have a value like: fi\ ~ J? m ~ 1, 
implying that it dominates the Universe just now! An easy way out would be 
to invoke anthropic arguments pj], but it is too feeble a justification, especially 
since the different contributions to A at the early phase transitions are 50 - 100 
orders of magnitude larger than what is observed. One would have to invoke 
a fine-tuning in order for the different contributions to cancel out, but yet not 
completely! However, theoretical models are being developed in an attempt to 
alleviate such problems (cf. |57[), some by allowing A to be a function of time 



(sec [148 and references therein) 
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